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The width of an x-ray energy level depends to a large extent on the probabilities of Auger 
transitions from that level. Consequently, rapid variations in line width with atomic number are 
to be expected whenever sudden changes in Auger transition probabilities occur. Atomic number 
ranges wherein large variations might be expected have been studied. Part of the width variation 
of L8; for elements 73=Z=92 is attributed to the Auger transitions Min—MyNiv,y. The 
transitions Ny—NinOu,101, NuNtv,vNiv,v, and Nin Nv, vNrv,y are shown to explain an 
intensity anomaly of M¢ for elements 38—=Z=58. Several other irregular line width variations 
are correlated with changes in the probabilities of Auger transitions. The evidence supports the 
theory that Auger transitions are primarily responsible for anomalous changes in line width. 





I. INTRODUCTION 


N 1935 Coster and Kronig! suggested that the 
relative diffuseness of the lines arising from 
transitions from L, initial states (hereafter called 
L; lines) was due to the broadening of the Ly 
level by Auger transitions. Large probabilities of 
Auger transitions from a given state decrease the 
mean life of atoms in this state and lead to an 
increase in the width of the energy level.? Theo- 
retical calculations for gold (79) by L. Pincherle* 
and Ramberg and Richtmyer* have shown that 
Auger transitions account for the greater part of 
the widths of all levels except K. The theory that 
the level width depended on Auger transition 
probabilities was confirmed experimentally by 
the work of the author,’ which showed that the 


1D. Coster and R. DeL. Kronig, Physica 2, 14 (1935). 
(1930) Weisskopf and E. Wigner, Zeits. f. Physik 63, 54 
*L. Pincherle, Nuovo Cimento (N.S.) 81, 162 (1935). 

*E. G. Ramberg and F. K. Richtmyer, Phys. Rev. 51, 
913 (1937). 
*J. N. Cooper, Phys. Rev. 61, 234 (1942). 


width of the Z; level increased rapidly with 
atomic number in the range 73<Z<81 where 
the probabilities for the Auger transitions 
Ly-Lin My and L;-LymMy are known to be 
swiftly increasing. 

Since there are many Auger transition proba- 
bilities which vary rapidly with atomic number, 
it is clear that there should be other levels (in 
addition to the Z;) which show sudden changes 
in width from one element to the next. It is the 
purpose of this paper to investigate the regions 
in which important variations in Auger transition 
probabilities occur and to correlate with these 
variations certain apparently erratic changes in 
level widths. 


II. PROBABLE AUGER TRANSITIONS 


According to Coster and Kronig,' the proba- 
bility of an Auger transition in which electron A 
goes from a state with a wave function ¥;(A) to 
a state represented by ¥2(A), while electron B 
passes from a state represented by ¥:(B) to one 
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TABLE I. Possible Auger transitions and the elements in which they may occur. 


























Lr>LuN; below Z=70 Mi—~Min Niv above Z=65 Ni Nin Ny1_ below Z=9 = 
Nu below 75 below 48 Nv below 53 Li 
Nin below 81 Ny above 65 Nivy below = 53 po 
Niv_ below 92 below 48 L 
NivNy _ below 84 L 
Lin M1 below 31 MivN1 below 88 Niv_ below 80 L 
Mu below 36 My below 35 L 
Mi below 37 My below 35 NvNy below 87 Li 
Miy above 73 L- 
below 50 LE 
My above 77 Min~MiwNy below 85 Num Nu1Our above 66 
below 50 Nrv_ below 84 below = 55 = 
Nin below 71 NiO below —s_& 
Nu below 65 above 66 
LumLin Mi below 30 N: below 55 below 55 
My below 30 My below 36 Nyu below = & = 
My above 90 My below 36 Nyt below g& | 
MyNy_ below 91 2s 
Niv_ below 89 NivNy | below 60 a) 
My MiuNvu below 91 Nin below 77 Niy below 60 IB. 
Nv: below 91 Nu_ below 70 L 
Ny below 72 Ni below 57 NyNy below © 
Ni below 71 My below 37 v 
Niu below 53 : ; Lp. 
Nu below 53 Miy>MyOu1 below 85 Ny NryvOr below 87 Lp 
J: below 47 Ou below 82 Nyur_ below &% . 
My below 34 Nyt below 85 = 
Mwy below 34 Ny below 57 
M-~Nun0m below 88 Nw below 5 thi 
MinNy below 87 Ou below 85 ‘ 
Ni below 50 O; below 79 NyOr below’ & the 
My below 34 Nvyu below 80 Nyu below 87 it i 
My below 34 Ny: below 80 Nv1_ below 86 
Ny below 53 Ny below 357 one 
MivMy below 44 Nrv_ below 53 Niv—>NyOw,v below — 8 bel 
Mwy below 43 ; 
elen 
MyMy below 45 Ny—>NvyiNyu_ below 9} T 
Ny: below 9% the 
NvuNvir below 9 calc 
— initi 
described by ¥2(B), is proportional to the square of W increases to a maximum and then falls off, — 
of the matrix element so that the probability of an Auger transition in z } 
which the ejected electron has a high energy is eB 
f ¥1*(A)W2(A)dripi*(B)p2(B)dre onal Zwei 
" rie Many Auger transitions are allowed for some 7 , 
elements and forbidden for others. An Auger r 
where dr; and dr2 are volume elements separated transition is possible only when the energy of a . 
by riz and the integration extends over all space. _ the initial state is greater than that of the final 7 
In order to obtain a large value of v, ¥:*(A) and _ state (of double ionization) and is most probable mek 
¥2(A) must overlap considerably as must ¥;*(B) when the energy difference is not too great. a 
and ¥2(B), while 7:2 should be small in the region Consequently, rapid fluctuations in the proba cal 
of appreciable overlapping. ¥i*(A) and y2(A) bility of such a transition are most likely in the liens 
generally overlap most when both initial and atomic number range wherein the transition first a 
final states of electron A lie in the same shell. becomes possible. In view of this fact it is d nenpns 
The energy W of the ejected electron is very considerable value to know the atomic numberat y 
important in determining the overlapping of which various Auger transitions first become per- 
¥i*(B) and ¥.2(B). If W is either very small or mitted. A systematic investigation was made of 
very large, the overlapping is slight and v will all Auger transitions which are allowed for some Fre 
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elements and forbidden for others. The result o 
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TABLE Il. Widths of L series lines in electron volts. 





























Line Transition Ag(47)* — Ta(73)— W(74)_—-RM(7S)-—Os(76)~—Ir(77)-—PA(78) = Au(79)-— TH(81)* ~—P(82)* U2) 
ip, LiMn 6.6 11.7 123 136 147 149 17.2 184 19.7 183 188 
IB. Mu 5.9 14.1 151 164 174 17.7 194 20.7 216 222 322 
= ie 11.0 11.6 121 132 140 146 — 181 193 199 394 
os Nu «10.2 10.0 110 120 126 133 156 164 178 194 324 
1B: Lu Mtv 2.4 6.6 6.5 6.5 6.6 6.8 8.3 8.5 9.0 84 14.3 
Ly. Nr 3.95 9.5 9.3 9.5 9.4 96 114 15 117 «49114 #159 
Ip. LuNv 3.72 9.8 9.8 9.8 9.9 9.7 i141 112 115 118 16.1 














«See reference 6 of text. % See reference 5 of text. ¢ See reference 9 of text. 4 See reference 7 of text. 


TABLE III. Ratios of line widths. 














Ratios Ag(47) Ta(73) Ww(74) Rh(75) Os(76) Ir(77) Pt(78) Au(79) TK81) Pb(82) U(92) 
LB3/LB: 1.77 1.20 1.26 1.39 1.49 1.54 1.59 1.72 1.80 1.55 1.17 
LB,/Lp2 1.59 1.44 1.54 1.67 1.73 1.83 1.82 1.96 2.00 1.88 2.00 
Ly2/LB2 2.96 1.18 1.24 1.35 1.41 1.51 — 1.66 1.71 1.69 2.45 
Ly3/LB: 2.74 1.02 1.12 1.23 1.27 1,37 1.40 1.49 1.58 1.64 2.02 
LB2/LB: 1.55 1.49 1.51 1.51 1.50 1.43 1.40 1.35 1.33 1.40 1.12 
LB:/Ly: 0.94 1.03 1.05 1.03 1.05 1.01 1.00 0.99 1.02 1.04 1.01 











this work is summarized in Table I, which lists 
the transition and the atomic numbers for which 
it is allowed. It will be observed that certain 
ones are allowed above one atomic number and 
below another, while they are forbidden for 
elements in between. 

The atomic number of the element for which 
the transition is first allowed was obtained by 
calculating the energy difference between the 
initial and final states. (A transition is obviously 
forbidden if the difference is negative.) Values of 
v/R for the various states were taken from 
Siegbahn’s Spectroscopie der Rontgenstrahlen, 
Zweite Auflage. The assumption was made that 
the energy of the final state of double ionization 
of element Z was equal to the sum of the energy 
of the inner level for element Z and the energy 
of the outer level for element Z+1. The energy 
so obtained is a good approximation, although 
not exactly correct. The uncertainty in the 
energy of the state of double ionization leads to 
an uncertainty in the atomic number at which a 
given Auger transition becomes possible. In the 
worst cases this atomic number may be incorrect 
by as much as three or four. 


Ill. WIDTH VARIATIONS AND 
AUGER TRANSITIONS 


From Table | it is clear that there exist many 
Auger transitions which are allowed for some 


elements and forbidden for others. Certain of 
these transitions unquestionably play an im- 
portant role in determining the widths of x-ray 
energy levels, while others are doubtlessly negli- 
gible. It is not immediately obvious which of the 
transitions become highly probable, although it 
would be possible to find out either by lengthy 
calculations or by detailed experiment. That 
there already exist sufficient line width data to 
establish the importance of certain of these 
transitions will be shown in this paper. 

There have been fairly complete studies of the 
widths of L series lines for several different 
elements. An exhaustive investigation of the 
entire ZL series of silver (47) was made by 
Parratt.* Measurements of the widths of the Z 
lines of uranium (92) are reported by Williams.’ 
The gold LZ series was measured by Richtmyer, 
Barnes, and Ramberg,* while Shrader® made a 
similar study of lead (82). Widths of selected L 
series lines for elements 73<Z<81 have been 
measured by the author.’ These data reveal 
many irregular changes in the widths of levels 
with atomic number. Most of these variations 
are readily explained in terms of Auger transi- 
tions listed in Table I. 


*L. G. Parratt, Phys. Rev. 54, 99 (1938). 
ok H. Williams, Phys. Rev. 37, 1431 (1931). 
* Richtmyer, Barnes, and Ramberg, Phys. Rev. 46, 843 


(1934). 
* R. E. Shrader, Cornell thesis (1936). 
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Fic. 1. Ratios of the widths of Z; lines to the 
width of L82(Zi11>Ny). 


In studying sudden variations in the width of 
a line as a function of atomic number one must 
remember that, in general, the width (in electron 
volts) increases with atomic number. This fact 
becomes of particular importance when one is 
considering elements that differ appreciably in 
atomic number. For example, the L series lines 
of gold (79) are all wider than the corresponding 
lines for silver (47), in most cases by a factor of 
two or three. This width increase for most lines 
is associated with the normal change with atomic 
number. One way of treating data which elimi- 
nates this problem is to deal with width ratios 
rather than widths. This essentially amounts to 
assuming that the ‘‘normal” width increase with 
atomic number is such as to leave the ratios 
unchanged. Available data indicate that this is a 
good approximation. 

Width ratios have one other great advantage 
when one is dealing with data taken by several 
observers, each of whom used different crystals 
on different spectrometers and then corrected for 
crystal diffraction patterns on fundamentally 
different assumptions. That advantage is that 
the ratios are far less sensitive to crystals, spec- 
trometers, and corrections than are the actual 
widths themselves. To be sure, the ratios are not 
independent of corrections for the finite resolving 
power of the spectrometers, but the percentage 
variation in the ratios is much less than that in 
the widths. 

In Table II are given the corrected widths in 
electron volts of seven L series lines for each of 
eleven elements, and in Table III the ratio of 
the width of each of these lines to the width of 
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LB: is recorded. Lf: is chosen as the referenge 
line because it arises from an Ly transition and 
the width of the Lim level is probably the least 
influenced by Auger transitions. The data op 
silver are by Parratt,® those on lead by Shrader! 
and those on uranium come from Williams’? 
paper. The line widths of the other elements were 
measured by the author.® It should be mentioned 
that the type of corrections applied by each of 
these observers was different. Williams’ data 
were corrected on the assumption that shapes 
were given by the Gaussian error function and 
his corrections were much smaller than thoge 
applied by the other observers. 

In Fig. 1 the ratios of the widths of LAs, Lg, 
Ly2, and Ly; to LB: are plotted as a function of 
atomic number. The first four lines all arise 
from transitions from the Z; state, while L, 
comes from an Ly—WNy transition. There js 
every reason to believe that the width of Lg, 
is behaving “normally” in this region; that is, 
neither the initial nor the final states are changing 
much in width. If this is true, it is clear that the 
width of the LZ; level must be increasing between 
tantalum (73) and lead (82) since all the lines 
arising from JZ; transitions show large and 
roughly equal increases in width in this range of 
atomic numbers. These width increases have been 
shown® to be associated with the rapidly in- 
creasing probabilities of the Auger transitions 
LyoLm My (possible above Z= 73) and L 
—LimMiv (possible above Z= 77). 

Since all the lines (except L8:) have the same 
initial state, the four curves should behave 
similarly unless the width of the final state is 
varying irregularly with atomic number. The 
curves do not behave alike above Z=81, so the 
conclusion is that here the widths of the final 
states must also be varying. In several cases the 
changes in width ratios are so great that the 
importance of certain Auger transitions are well 
established, while in other cases the variations 
are only of the order of the experimental errors. 
In such cases the variations may be due to 
Auger transitions and have been attributed to 
them, although it is possible that they are due to 
a fortuitous combination of possible errors. 

The line Lf; narrows greatly between thallium 
(81) and uranium (92). Since L; arises from aa 
I;—Mim transition, this means that the Mm 
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level must be decreasing in width. Table | re- 
yeals that the Auger transitions Mi—MyNiv 
and Mir—MyNy become possible for elements 
of atomic number below 90 and are forbidden 
for elements of higher atomic number. Also the 
transitions Mim—MivNiv and Mir—MiyNy are 
forbidden for uranium and allowed for the other 
elements. (See Fig. 2, reference 10 and Fig. 1, 
reference 11.) Judging by the great drop in the 
relative width of L§; these transitions must ac- 
count for much of the width of the My level for 
elements around Z=81. 

Since the probability of the transitions will 
decrease as the energy of the ejected electron 
becomes too great, the width of the My level 
might be expected to be less for Ta(73) than for 
T1(81). This argument is supported by the obser- 
vation that the actual width increase of LA; 
between these elements is greater than that of 
any of the other L; lines. This would be expected 
if the Auger transitions Mm—MyNjy and 
Mm—MyNy had a maximum probability in the 
neighborhood of Z=81. 

Further support for this explanation can be 
found in the work on Ma-satellites by Hirsh,’° 
who has shown that these satellities have a 
maximum relative intensity near lead (82). Since 
these satellites arise from an MyNphy.y initial 
state, they will obviously increase in intensity as 
the Auger transitions Min—MyNjy,y become 
more probable. Thus the importance of these 
radiationless transitions is established by inde- 
pendent observations. 

The width of L8,(Z:—>Mu) appears to be rela- 
tively stable between thallium and uranium al- 
though data on intermediate elements might 
indicate width variations associated with the 
Auger transition My—MjyN; which is possible" 
for elements of atomic number below 88. 

Both Ly2(Li1—Nn) and Ly3(Z:1—Nim) are much 
wider for uranium than might have been ex- 
pected, showing that the Nn, and particularly 
the Nin, levels have become very wide. There is 
nothing in Table I which would predict such an 
increase. The addition of ten electrons in outer 
levels between lead and uranium will make many 
new Auger transitions possible and it is con- 
ceivable that one of them might be very probable. 


’F,R, Hirsh, Phys. Rev. 57, 662 (1940). 
4 F.R. Hirsh, Phys. Rev. 59, 766 (1941), Fig. 1. 
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Figure 2 shows the width ratios L6:/L8,; and 
LB,/L7:. LB: arises from the transition Ly;— Ny, 
while Lf, and Ly; have Ly as an initial state and 
their final states are Myy and Nyy, respectively. 
It is immediately clear that these ratios are much 
more constant than those involving LZ; lines. The 
width ratio LB2/Ly,; is almost constant for all 
the elements listed, but the ratio L8:/L8; de- 
creases with atomic number. There are a number 
of Auger transitions which may play a part in 
this. First, 18, may be wider at uranium because 
the Auger transition Ly—Zi1My becomes pos- 
sible above Z=90 and may contribute to the 
width of the Zy level. If this is true, the width 
of Ly: should be similarly influenced, and in 
order to explain its constancy, it is possible to 
point out that the Nyy level may have narrowed 
because the transitions Nry—NyOry.y are for- 
bidden for elements above 81. Part of the decrease 
for the width ratio L§,/L8,; may be due to a 
narrowing of the Ny level for uranium, since the 
transitions Ny Ny Nyivu and Ny Nyu Noi 
are impossible for that element. Another factor 
which contributes to the small ratio for uranium 
is the smaller correction subtracted from the 
widths in this case. (Obviously, too small a 
correction would tend to make this ratio nearer 
to unity.) 

Some interesting changes in the widths of 
levels occur between elements 73 and 47. For 
silver the lines Ly: and Ly; are nearly twice as 
wide as L8; and Lf, which indicates that the 
Nu and Nyy levels are much wider than the My 
and My levels. For elements 73 to 82 the My 














iso @—~»—?— 
> 
. lied aie 
= on da 
ae a a re” lee ie Pee 
- OOOO gm nnn menmnnn mn nn 4 
1) 
Sy 
S 
3 

0.50 . 
s 
2 ° yg © Ne 

° i 4 i rn 
7) 74 78 82 86 90 
ATOMIC NUMBER 


Fic. 2. The width ratios L8:/L8; and Lf;/Ly,; as a 
function of atomic number. (Transitions: L8:, Lin—~Ny; 
Lpi, Lu My; Lyi, Lu Nv.) 
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TABLE IV. Widths of M¢ in electron volts (after Kiessig). 











Element Width Element Width Element Width 
Sr(38) 1.99 Rh(45) 4.57 Sb(51) 18.5 
Zr(40) 2.06 Pd(46) 6.74 Te(52) 18.5 

Nb(41) 3.44 Agi?) 13.6 Ba(56) 7.0 

Mo(42) 3.05 Cd(48) 14.8 La(57) 9.1 
Ru(44) 3.79 Sn(50) 16.4 Ce(58) 10.6 








and My levels are wider than the corresponding 
N levels, but for uranium the WN levels are again 
wider. This can be partially explained as follows: 
(1) For elements of atomic number 70 to 82 the 
Mi level should be wide because of the high 
probabilities of the transitions Miy—MyMiv,v 
and (possibly) Mir—MtyMiv,v which have al- 
ready been mentioned in connection with the 
small width of UL. Also, for silver the Auger 
transitions Mj;—MiyN; and Mm—MyNji,n are 
impossible, whereas they should contribute to 
the width of the Mm level for elements of higher 
atomic number. (2) The Mn level may be ex- 
pected to be wide for elements not too far above 
70 because of the Auger transitions My— Min Mv 
and My—MinNy. These transitions are for- 
bidden for elements between (about) 65 and 48. 
Whether or not they are allowed for silver is 
doubtful, but judging by the small width of the 
Ag Mn level it seems certain that these transi- 
tions are very improbable for silver even if they 
are possible. 

While these Auger transitions may serve to 
explain why the M levels are wider than the N 
levels for elements of atomic number between 73 
and 82, still the widths of the Ny and Nm levels 
for silver are so very great that their excessive 
widths require further explanation. The Ny level 
is widened by the transitions Nno—-NmOun,m 
which are possible for elements of atomic number 
below 55 and by Nu-—Mrz,vMv.v which are 
allowed below Z=60. The great width of the 
Nin level is attributed to the transitions Nin— 
Nrv,vNrv,v which are possible for elements below 
Z=56. 

Line width data by Kiessig’? support the con- 
tention that these transitions must be highly 
probable for silver (47). In Table IV are the 
widths of M¢ for fifteen elements as reported by 
— The line Mg arises from a transition 


Kiessig, Zeits. f. Physik 109, 671 (1938), and 95, 
555 loge 
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My—Nu,m. The width increases rapidly from 
Z=38 to Z=52, showing that the Ny and Min 
levels are widening rapidly. (The My level js 
“normal” for silver and there is no reason to 
think that these width changes are due to the 
initial level.) The width has decreased tremen- 
dously at Z=56, just as would be expected from 
a consideration of Auger probabilities. Kiessig 
did not offer an explanation for the intensity 
anomaly, but pointed out that the elements at 
the beginning of the fifth period of the periodic 
table (36-54) had small M¢ widths and those at 
the end had great widths, while the three ele. 
ments at the beginning of the sixth period again 
had narrow breadths. Siegbahn and Magnusson” 
had previously found that the Mf¢ lines of Br(35) 
and Rb(37) were very sharp, so that the general- 
ization that the width increased as one went 
through an atomic period was not justified. Now 
it seems virtually certain that the above-men- 
tioned Auger transitions provide the correct ex- 
planation for this intensity anomaly. 

It is not inconceivable that the width increase 
between element 56 and 58 is due to the in- 
creasing probability of -Ni—Niv,vNiv.v but the 
data do not cover a big enough atomic number 
range to be conclusive. If the width of M¢ were 
shown to be much less above Z = 60, this explana- 
tion would be much better established. 

The width of the Ny level is consistently wider 
than that of the Ni level for all elements con- 
sidered here. This is in accord with the general 
rule that for different levels having the same 
total quantum number the level width decreases 
with increasing azimuthal quantum number. An 
exception to this rule is found in silver where the 
Myx level has a greater width than the My level. 
The unexpectedly large Min width may be associ- 
ated with relatively great probabilities of the 
transitions My—MiyN; and Mi—MyN,, pos- 
sible below atomic numbers 55 and 57, respec- 
tively. 

It is interesting to note that the width ratios 
LB2/LB8; and LB2/Ly; for silver are almost the 
same as for elements of atomic number 73 to 81. 
Since all the levels involved in these lines are 
associated with relatively high azimuthal quan- 
tum numbers, the Auger transitions generally 


13M. Siegbahn and T. Magnusson, Zeits. f. Physik 8, 
559 (1934). 
































where 
from t 
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contribute much less to the level widths. The 
relative constancy of the width ratios for these 
lines over SO great an atomic number range sup- 
ports the ratio method for handling width data 
for elements of vastly different atomic number. 


IV. CONCLUSIONS 


Since the width of a state is inversely propor- 
tional to the mean life of atoms in that state, a 
very probable Auger transition increases the 
level width. This results in a width increase for 
all lines which have this level as either an initial 
or final state. Because Auger transition proba- 
bilities frequently vary rapidly as a function of 
atomic number, one should expect rapid changes 
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in line width. This prediction is verified. In this 
paper it has been shown that the width variation 
of L8; for elements 73<Z<92 can be accounted 
for in terms of the Auger transitions Miy— 
MyNry,y. Also the anomalous width behavior of 
Mf reported by Kiessig has been shown to be 
associated with the rapid changes in the Auger 
probabilities No—NmOnm, Nu—-Ni,vNrw,v, 
and Nim Nv, v Nrv,v for elements 38<Z<58. 
Certain other irregular width variations have 
been shown to be explicable in terms of Auger 
transitions. More data on line widths would 
doubtless reveal other cases in which line widths 
are increased by Auger transitions. 

The author wishes to thank Dr. F. R. Hirsh, 
Jr., for his encouragement and suggestions. 
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On Apertures of Transmission Type Electron Microscopes Using Magnetic Lenses 


L. Marton AND R. G. E. Hutter 
Division of Electron Optics, Stanford University, Stanford University, California 
(Received January 5, 1944) 


The numerical aperture of light microscopes is larger than that of electron microscopes at 
their present state of development by a factor of 100 to 1000. This is due to the fact that 
electron lenses have higher aberrations than the highly corrected glass lens systems used in 
modern light microscopes. In the design of an electron microscope care has to be taken therefore 
that the proper numerical aperture should be used to give the optimum resolution. The majority 
of transmission type electron microscopes use at present magrietic lenses. For magnetic lenses 
of the axial field distribution H(z) = H»/[1+(z/a)*] optimum conditions for size and location 
of apertures will be stated. The behavior of the condenser lens-objective lens system with 
respect to the angular aperture of the illuminating electron beam also will be discussed. 


I. THE MAGNETIC LENS OF THE AXIAL FIELD 
DISTRIBUTION H(z) = H/(1+(z/a)?] 


HE small values of numerical apertures of 

electron microscopes can be achieved by 
allowing only such electrons to take part in the 
image formation whose paths do not pass through 
the lenses too far away from the optical axis. 
These paths are called paraxial. They can be 
described mathematically as the solutions r = r(z) 
of the so-called paraxial ray differential equation: 


d*r(z) 
dz* 





é 
2 - 1 
oe al (z)r(z) =0, (1) 


where r= r(z) signifies the distance of the electron 
from the axis at the point z—the axial coordinate; 


H(z) is the axial field distribution; e/m is the 
specific charge of the electron, and V is the 
accelerating potential for the electrons. 

It has been shown! that for an axial field dis- 
tribution, very closely approximated in practical 
magnetic lenses, 


H(z) = Ho/(1+ (2/a)?). (2) 
The solution of Eq. (1) is given by 
r(z) =a(1+(z/a)*)! 
X {Ci sin [(1+K?)! are cot (z/a) ] 
+(C: cos [(1+K*)' arc cot (z/a)]}, (3) 


1 Walter Glaser, Zeits. f. Physik 117, 285 (1941). 
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where Ho is the maximum field strength on the 
axis, z=-+ta is that coordinate where H(z) 
= H(+a)=}H, (2a is termed the half-width of 
the field curve), and 


K?=eH,a2/8mV (4) 


is a parameter characterizing the lens strength. 
C, and C; are two arbitrary constants of integra- 
tion, which can be determined so that the path 
satisfies any initial or boundary conditions. 


II. ANGULAR APERTURE. THE CONDENSER 
LENS-OBJECTIVE LENS SYSTEM 


The numerical aperture is defined by 
N=n-:sin a, 


where m is the index of refraction and a is the 
so-called angular aperture, that is, the largest 
angle formed by light rays which pass the 
center of the object. In electron microscopes 
n=1 and sin a~a for the values a=10~ to 107 
radian. 

Borries and Ruska? calculated the angular 
aperture of a lens system, consisting of a con- 


Condenser lens 
aperture 







=—— Object 
a al 











4nin a cee 
Minimum cross-section 


Fic. 1a. Condenser lens. Notation for Eq. (5). 

















Fic. 1b. Angular aperture in function of the refracting 
7 a the condenser lens according to B. v. Borries and 
. Ruska. 


2B. v. Borries and E. Ruska, Zeits. f. tech. Physik 20, 
225 (1939). 
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denser lens and an objective lens, for various 
lens strength of the condenser lens (Fig. 1a), 
Their results are given by the following formulas 
and graphs (Fig. 1b). 


dmin/2 
a= 
a+b— (a-b/f.) 



















for 0<-< 


1 a+b dmin 
fe \a-b gS 





a 
Ac = Ac max >= — A, 


a+b dmin 1 a+b dmin 
for ({——-—- )<-<( + ), 






(5) 














a:b 2a’a, a-b 2a?’a, 
dmin/2 
a, = 
(a-b/f.) —(a+6) 















for 


a+b Gate 1 
— ) <-<e, 

a-b 2aa/ f. 

where dmin is the diameter of the minimum 
cross section of the electron beam in front of the 
cathode (the virtual source), a is the distance of 
this minimum from the condenser lens, and } js 
the distance of the object from the same lens, 
f. is the focal length of this lens and a, is the 
angle given in Fig. la. 

In these calculations, of purely light optical 
character, it was assumed that no other lens was 
located between the condenser lens and the 
object. This, however, is no longer true if the 
object is located so that it is well within the field 
of the objective lens. Glaser! has shown that for 
the condition of minimum spherical aberration 
this location is imperative. In this case a part of 
the field of the objective lens lies between the 
object and the condenser lens and has a focusing 
effect on the incoming illuminating electron 
beam, changing the angular aperture as given by 
the condenser lens strength alone. 

Assuming for both lenses a field form given by 
Eq. (2) it is possible to recalculate the value of 
the angular aperture, taking into account the 
focusing action of the objective lens fragment. 
It is necessary, however, to assume a field free 
plane midway between the objective lens and 
the condenser lens—an assumption normally 
realized in electron microscopes (Fig. 2). 






























the | 
objec 
the ¢ 
pass 


APERTURES 


OF TRANSMISSION TYPE 


ELECTRON MICROSCOPES 163 


Field Free 
midplane 


' 
' 
‘ 
' 
' 
' 
' 
' 
i 
4 





Condenser lens 


Ha) . ates 


i+(ah 


Notation: 
Hoo 


2a 
2 
2 eHoade 
oo = 2 
8mV opal 
or 2 
p we = (1+ Ko)? 


A X4e-D 
Objective lens 


.* 
4 (zy 


Hiay« 


He } maximum field strengths of the lenses. 


goo half-widths of the field strength curves. 


parameters characterizing the strengths of 
the lenses. 


Fic. 2. Electron path through the condenser lens and objective lens. 


The values of w realized in practice are 
1=w=2;r, is the radius of the cross-over disk of 
electrons in front of the cathode of the electron 
microscope; R, is the radius of the condenser lens 
aperture; R, Ri, r2 are ordinates of the electron 
path at the abscissas z= —D/2, 2=0, z=+a, 
respectively. 

The method for arriving at a solution for the 
angular aperture a consists of the following steps. 
Two solutions of the paraxial ray equation, one 
for each lens, are taken and matched at the field 
free midplane so that their distances from the 
optical axis are equal and their slopes are iden- 
tical. The total path has to pass the center of the 





object and must also go through a point on the 
circular disk of the electron cross-over. Among 
all such paths the one is taken which gives the 
largest angle a for given lens strengths and 
position coordinates. This amax is then the angle 
determining the angular aperture. 

Glaser! has shown that for minimum spherical 
aberration the object would have to be very near 
the center of the objective lens. For the sake of 
simplicity it is assumed here that the object is 
at the center of the lens. 

Rewriting the formulas of Borries and Ruska 
in terms of quantities used here, we find the 
angle a@max is given by 


—TaWe 





tan Qmax 


By varying w,. between 1 and 2 a variation of 
@max is obtained which is identical with that given 
by the previous formulas. 

If w. becomes 


T 





O.= ’ (7) 
arc cot (a/a.) —x+arc cot (D/a.) 

the cross-over would have its image at the 

object position, which means that rays leaving 

the center of the cross-over at any angle would 

pass through the center of the object at any 


~ ae{[1+(a/a.)*1[1+(D/a.)?]}! sin {w.Lare cot (a/a.) —x-Farc cot (D/a.)]} 


(6) 





angle up to 90°. This angle, however, cannot 
be reached because of the limiting effect of the 
physical aperture of the condenser lens (Fig. 3). 


Condenser lens 





Minimum cross-section 


( radias = r. ) 


Fic. 3. Limiting effect of the physical aperture of the con- 
denser lens. 
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There will be a value we: which realizes path (1) and a greater value we2(we2>we:) which 
realizes path (2). In order to find the limiting w, it is necessary to find the distance Ri, at which 
any path passes the condenser lens center (where the physical aperture is assumed to be located), 


R, is given by 
a¢ Dy *7 Dr 
Ri=“]1+(=) | tan a sin [(are cot —-)| (8) 
We ae ae 2 
together with 


Qe a\? Dv21)3 a D 
r.=— [1+(=) +) |} tan a@ sin [o(»-are cot ——arc cot —)| (9) 
We a- ae ae ae 


the ratio R:/r, can be plotted as a function of w,: 


Ri 1 sin {w.[arc cot (D/a.) —(x/2) ]}} 
tr. e [1+(a/a.)?}! sin {w.[*—arc cot (a/a.) —arc cot (D/a.)]} (10 













we, and weg are the values for which Ri/r,=+R,/r-. For all values of we, we: Sw,Swee the physical 


aperture becomes the limiting factor for the angular aperture. 
Repeating now the same procedure for the case where the fragment of the objective lens is active 


in forming the angular aperture, we obtain the following expressions: 
(1/a.)rewel 1+(D/2a.)*}*/[1 + (a/a.)*}* 
m (K1/a.)w. cos {w.[#*—arc cot (a/a.) —arc cot (D/2a.) ]} 
+ {(K:D/2a.*)+K2[1+(D/2a.)*]} sin {w.[*—arc cot (a/a.) —arc cot (D/2a.)}} 






» (If) 








tan Qmax 















where 
a D\*} rg D 
K.=“]1+( ) sin [eo(=—are cot — | 
w 2a. 2 2a. 
1 1 D , D ba D 
K2=— — sin [eo( =—arc cot — +. cos [oo( =—arc cot —)|}. 
wo [1+(D/2a,)?}*\ 2ao 2 2a. 2 2a 






K, and Ky depend on wo, the parameter of the objective lens, and the location of the object; these 


quantities are plotted in Fig. 4. 
In order to determine the range of w, values for which the physical aperture of the condenser lens 


becomes the limiting factor for the angular aperture a the distance R, must be determined. 
The equation giving R, in this case is 


a, 1 We T D 
R,=— K, = os | we( =—are cot )] 
we [1+(D/2a,.)*}* la, 2 2a. 


D &K; D\? 7 D 
+|—+=(1+(—) )| sin [e(=—are cot —)|I, (13) 
2a,? K, 2a. 2 2a. 
and from Eq. (11) 


af 1+(a/a.)? PP [we a D 
of Kii— cos [+.( x—are cot ——arc cot — 
wel 1+(D/2a,)? ae ae 2a. 


D Kz D\? a D 
+ +=(1+(—) ) ]sin [.(x-arc cot ——arc cot )]| tana. (14) 
2a.2 Ky 2a. a. 2a. 
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Forming again R;/r, and plotting this ratio as 

a function of w, we find the two values «,, for 
which Ri/re= +R-/te. These are called wes, wes. 
For all values w for which w.Sw.Sw.2 the 
angular aperture is determined by the physical 
condenser lens aperture. 

Furthermore, there will be one value w=w* 
for which the denominator in Eq. (11) -becomes 
zero (or Ri/re=+*). For this value w* the 
cross-over of the electrons in front of the cathode 
is focused on the object. This means that all rays 
starting at the center of the cross-over will pass 
through the center of the object. In Eq. (11) ra 
and the denominator become zero; hence tan a 
is indeterminate, as it should be. 

Figures 5 and 6 are obtained if the following 
values are assumed: 

D=420 mm, a= 250 mm, a,=3 mm, do9=1.5 mm; 
r-=0.05 mm, R,=0.4 mm, hence R./r,=8; 
wo= 1.5, isa. 





| l J | 
1 1.2 1.4 1.6 1.8 


Fic. 4. K; and K¢ in function of w. 
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Fic. 5. Angular aperture of the incoming beam in function 
of the refracting power of the coridenser lens. 































FiG. 6. Ri/re as a function of we. [wer = 1.0059; w* = 1.0061; 
wee = 1.0063. ] 
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Fic. 7. Location of the aperture on the optical axis in 
on of the magnification and the lens parameters a 
and w. 


For the sake of comparison, Fig. 5 shows amax 
as calculated by means of the indicated methods 
and as calculated for the same numerical values 
by the method of v. Borries and Ruska (dotted 
line). 

Contrary to previous assumptions it is shown 
that the angular aperture cannot be decreased 
below any prescribed limit by increasing the 
refracting power of the condenser lens. There 
exists a definite minimum for a, as observed in 
electron microscopes, because of the effect of 
refracting properties of the objective lens 
fragment. 

The formulas derived so far allow the deter- 
mination of the angular aperture of the incoming 
electron beam. After passing the specimen the 
electrons are scattered and diffracted causing 
the angular aperture to be increased. In order to 
limit the angular aperture to values prescribed 
by limitations of the resolving power physical 
apertures must be used. These apertures are 
usually located at the center of the objective 
lens. It is obvious, however, that the diameter 
of the aperture must become smaller and smaller 
the nearer the object is to the center of the lens, 
if it has to achieve the same ray limiting action. 

Knowing the exact solution for the electron 
path it is possible to determine an optimum 
position and size for such an aperture which 





AND R. G. E. HUTTER 






achieves the same ray limiting effect. Let 
y(z)=r(z)/a and x=2/a; (15) 
then Eq. (3) can be written as 
y(z) = (1+x*)!{ C, sin [w arc cot x ] 
+C;cos[warccot x]}. (16) 


If the constants of integration C; and C; are 
determined so that the path goes through a 
point x=xo, the x coordinate of the object posi. 
tion, and has there a slope y’(xo) =tan a, Eq. (3) 
becomes: 


1 
y(x) =—[(1 +2?) (1+?) }! sin [w(arc cot xo 


—arc cot x)]tana (17) 


or 
sin [w(arc cot x»—arccotx)] tana 





y(x) = (18) 


sin [arc cot xo | sin [arc cot x] w 
The slope of this path at any point x is given by 
tan 5 al 
w L 1+x? 


X {x sin [w(arc cot x»—arc cot x) } 





yy’ (x0) = 


+w cos [w(arc cot x»—arc cot x) ]}. (19) 


Let x, be the x coordinate of the point of 
maximum deviation of the electron from the 
axis. Then x, is given by the transcendental 
equation: 

° xr 1 Sen 
arc cot X»=arc cot X»+——— arc cot —. (20) 


®o ®W 7) 
The next step is to combine Eq. (8) with the 
formula for the magnification of the lens.' 


1 1 





. : (21) 
xo—cot (x/w) sin (x/w) 


M is then given as a function of x,, and the lens 
parameter 





_ sin {are cot %m+(x/w) —(1/w) arc cot (%m/w)} (22) 





For large values of M (M>100) Eq. (22) 
becomes: 


sin (x/w) 
M =3————_-x,,,*. (23) 
(1—w?) 


sin {(1/w) arc cot (%»/w) —arc cot Xm} 





The position coordinate of the maximum path 
deviation from the optical axis is plotted in Fig.7 
as a function of the magnification for several 
lens parameters. The value w= 1.06 corresponds 
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to a weak lens; (w= 1 characterizes a lens of zero 
refracting power!) and w=2 to a strong lens. (It 
has been shown by Glaser' that for values #=2 
an object has two images.) 

An aperture placed at x=x,, on the optical 
axis can have a considerably greater size than 
one placed at the center of lens and yet have the 
same ray limiting effect. The size of the aperture 
for any angular aperture a@ (as determined for a 
permissible spherical aberration) can be deter- 
mined by the use of the formula: 


radius of aperture fa 








tan @ 
1+xm?)(1+x 97) 7! 
4 +xm?)(1+x . (24) 
(w?+%m”) 


tan a 





2r,/a tan @ is plotted in Fig. 8 as a function of 
the magnification for various values of w. 

The ordinate of the path at the center of the 
lens is given by 


a 7 
r(0) =—(1+20?)! sin [o( arc cot »—")| tan a. 
w w (25) 


The ratio of 2r4/2r(0)=da/do is plotted in 
Fig. 9 as a function of the magnification for 





Fic, 8, Reduced aperture size in function of the magni- 
fication for various lens strength parameters. 
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Fic. 9. Ratio of aperture sizes of two equivalent aper- 
tures, one located at optimum position (diameter d,4) and 
the other located at the center of the lens (diameter do). 


various values of w. It shows the gain in 
aperture size for an aperture at the optimum 
position as compared with the customary aper- 
ture design at the center of the lens. The ratio is 
given by 





TA 1+,” 4 1 
ee cs 


r(0) w*+x,7J sin {w[Larc cot xo—(x/2)]} 


The stronger the lens the more unfavorable is 
the aperture position at the center of the lens. 

Because of the increased size, the difficulties in 
machining the apertures are reduced. Another 
advantage results from the location of the 
aperture outside the lens inasmuch as the 
removal for cleaning purposes is facilitated. 
Furthermore, the exact location of the aperture 
on the axis is not too critical since the path is 
parallel to the optical axis at the optimum 
position and hence changes its slope less rapidly 
than at other positions. 
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The Ionization Cross Section of the Silver Z:; State 


J. J. G. McCue* 
Cornell University, Ithaca, New York 


(Received December 31, 1943) 


The ionization cross section of the silver L111 state has been measured, in arbitrary units, by 
observing the intensity of the La doublet radiation from a thin target of silver bombarded by 
cathode rays whose energies ranged up to nine times the excitation energy of the Zit state. 
The observed data are corrected for the effects of diffusion, rediffusion, and retardation of the 
cathode rays. The measured values of the cross section agree fairly well with calculations based 
on the Born approximation, but the agreement is not within the estimated experimental error. 








HE dependence of the ionization cross sec- 
tion of an atom on the energy of the im- 
pinging electron is a matter that has received 
considerable attention since it affords a means of 
testing the applicability of wave mechanics to 
collision processes on the atomic scale. One can 
measure the ionization cross section of the inner 
shell of an atom by observing the intensity of an 
x-ray spectrum line whose initial state involves 
ionization of the inner shell in question, the 
ionization being produced by subjecting the atom 
to cathode-ray bombardment. The chief difficulty 
arises from the necessity of using monokinetic 
cathode rays, but if the target be made in the 
form of a film, then it can be made thin enough so 
that the cathode rays pass through it with only a 
small energy loss, for which a correction can be 
made. Using such targets, Webster'-* and his 
colleagues have measured the ionization cross 
section of the silver K shell for electrons whose 
energies ranged from one to seven times the bind- 
ing energy of the K electrons. On the publication 
by Burhop*‘ of calculated values of the ionization 
cross sections for the silver L states, it seemed 
worth while to measure the cross section for 
the silver Ly state in arbitrary units by using 
the thin target techniques originated by Webster. 























I. THE APPARATUS 


The x-ray tube designed especially for this 
investigation is described elsewhere.® 







* Now at the Radiation Laboratory, Massachusetts In- 
stitute of Technology, Cambridge, Massachusetts. 

1 Webster, Clark, Yeatman, and Hansen, Proc. Nat. 
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The tube voltage was obtained from a motor 
generator set coupled to a high voltage trans. 
former and full-wave kenotron rectifier, with a pj- 
section filter composed of two 0.25-microfarad 
condensers and a large choke. Tests with a 
cathode-ray oscilloscope failed to reveal any 
ripple in the high voltage under operating con- 
ditions although a ripple of 50 volts would have 
been detected. 

A Siegbahn-Thoraeus® vacuum spectrometer, 
with a selenite crystal and an ionization chamber, 
was coupled to the x-ray tube by means of a brass 
pipe whose axis cut the axis of rotation of the 
crystal and the center of the target, making an 
angle of 5°+1° with the plane of the target. The 
pipe, which had an inside diameter of 0.5 inch, 
could be blocked by a movable lead stop. The 
stop was moved into the beam from time to time 
during the intensity measurements to keep track 
of the zero reading of the galvanometer con- 
nected to the detecting circuit. Two lead jaws 
fixed on the front of the ionization chamber 
limited the width of the beam entering the 
chamber. The other limiting slit in the system 
was the focal spot itself. Viewed from the crystal, 
that is, at an angle of 5° from the plane of the 
target, the focal spot had a width of the order of 
0.5 mm and a height not exceeding 6 mm. 

The selenite crystal was 10 mm high and 25 
mm long, with its axis of rotation 95 cm from the 
focal spot and 15 cm from the slit at the ioniza- 
tion chamber. This slit, 14 mm high, was 1.65 mm 
wide, the width being chosen so as to allow a band 
of radiation 23 x.u. wide to enter the chamber 
when the spectrometer was set on the Ag La line 


*M. Siegbahn and R. Thoraeus, J. Opt. Soc. Am. 13, 
235 (1926). 
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at 4150 x.u. This ensured the collection, in the 
chamber, of all of the Ag La doublet radiation 
which the crystal reflected in the first order.’ 
Moving the ionization chamber through 3 min- 
utes of arc on either side of the peak of the line 
decreased the intensity reading by only 1 percent. 
Therefore slight changes in the position of the 
focal spot could not introduce any error into the 
measurements, for the focal spot would have had 
to move 4 mm to shift the peak of the line 3 
minutes. 

The ionization chamber that took the place of 
the plate holder with which the spectrometer was 
originally equipped was a hollow brass cylinder 
12 cm long and of 5 cm inside diameter, with a 
collecting rod mounted parallel to, but displaced 
from, the axis of the cylinder. The chamber con- 
tained air at atmospheric pressure. Checking the 
variation of ionization current with x-ray tube 
current, with the x-ray voltage held constant, 
showed that in the range of intensities encoun- 
tered, the ionization current was proportional to 
the x-ray intensity within the accuracy of the 
observations, which was about 1 percent. 

The ionization current was amplified by an 
FP-54 electrometer tube operating in a Barth® 
circuit connected to a Leeds and Northrup type R 
galvanometer whose sensitivity was 4xX10-'° 
amp./mm/meter. The over-all sensitivity of the 
amplifying system was about 30,000 mm/volt 
with the scale 170 cm from the galvanometer. The 
non-linearity of its response was detectable, but 
the range of intensities used was so small that the 
error introduced by assuming a linear response is 
less than 0.5 percent. During the course of the 
measurements, the sensitivity varied somewhat 
because of changes in the temperature of the 
battery. Frequent measurements of the sensi- 
tivity made possible appropriate correction of the 
data. 

The target consisted of a silver film distilled 
onto a beryllium disk,* 24 mm in diameter and 
0.8 mm thick, soldered with Alumaweld to a 
water-cooled beryllium backing. The interference 
pattern formed when a glass optical flat lay on 
the disk showed that the target surface was flat 


*L. G. Parratt, Phys. Rev. 54, 99 (1938). 

*G. Barth, Zeits. f. Physik 87, 399 (1934). 

* Beryllium disks were generously furnished by the 
Brush Beryllium Company. 
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within one wave-length of green light over the 
area covered by the focal spot, while the inter- 
ference rings showed that the smoothness of the 
beryllium surface was about the same as that of 
plate glass. 

The silver, 99.98 percent pure, was melted in a 
horizontal tungsten trough 10 cm beneath the 
horizontal polished surface of the beryllium disk ; 
an ionization gauge showed that during the dis- 
tillation the pressure in the chamber was 3 X 10-5 
mm of Hg. Beside the beryllium surface lay a 
clean, clear-glass microscope-slide. A movable 
baffle plate allowed the beryllium and glass sur- 
faces to be exposed simultaneously, for an appro- 
priate time, to the silver coming from the trough. 
Since the heating current in the trough was 
turned on for only about 2 seconds, the silver 
condensed on a surface whose temperature was 
near 20°C. 

The thickness of the film was obtained by 
measuring, at three wave-lengths in the visible 
spectrum, the opacity of the film on the micro- 
scope slide. The data published by Strong and 
Dibble’® make possible an estimate of the film 
thickness when the opacity at these waye-lengths 
is known. The surface density of the target film 
according to their curves was 23 micrograms per 
square centimeter, but the data they obtained by 
direct weighing indicate that the actual surface 
density is greater than that shown in their curves ; 
giving equal weights to their two methods of 
weighing leads to an estimated density of 27 
micrograms per square centimeter for the target 
film. It is sufficient to assume that the volume 
density of the film was the same as that of 
massive silver, which implies that the thickness 
of the target was 260A. A thinner film would have 
been desirable in some respects, but the thickness 
of the thinner film would be so uncertain that 
260A seemed to represent an optimum thickness. 

During the measurements, some contamina- 
tion of the silver by material from the filament 
occurred. The portions of the film that could 
“see’’ the filament became discolored, but the 
focal spot was not visibly affected. Tests on the 
excitation potential of the silver La lines showed 
that the cathode rays did not experience any 


%” J. Strong and B. Dibble, J. Opt. Soc. Am. 30, 431 
(1940). 
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Fic. 1. The observed radiation and the continuous 
ckground. 








serious retardation in a surface film before 
striking the silver. With a massive silver target 
that had been used for some time and was very 
much dirtier than the thin target ever became, 
La radiation could be detected when the tube 
voltage was only 150 volts above the Ag La exci- 
tation potential (3.35 kv). The surface film on the 
massive target therefore retarded 3.5 kilovolt 
electrons by not more than 150 volts; the re- 
tardation was of course less for cathode rays of 
higher energy. 

The safe limit of power. dissipation in the 
target was about 30 watts. In order to avoid 
possible injury to the polish of the beryllium 
surface, it was touched after the polishing opera- 
tion by nothing but lens paper. There seems to 
have been enough grease on the beryllium to 
inhibit good thermal contact between the silver 
and the backing. 


Il. THE DATA 


After tests had been performed to make certain 
that the deflection of the galvanometer was pro- 
portional to the tube current at constant voltage, 
and was independent of the size of the focal spot, 
the intensity of the La doublet (4.15A) emitted 
from the film was measured at various voltages 
from 4 to 30 kv. A set of auxiliary measurements 
made with the spectrometer set at 4.30A de- 
termined the intensity of the continuous radia- 
tion in that part of the spectrum. The results 
appear in the upper and lower curves, respec- 
tively, of Fig. 1. The ordinates are centimeters of 
galvanometer deflection per milliampere of target 
current, which will be called the “specific in- 
tensity.”’ The scattering of the points at the right 
end of the upper curve is due to unsteadiness in 
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Fic. 2. Specific intensity of the La doublet. I—uncor. 
rected; II—corrected for fluorescence; 11I1—completely 
corrected. 





the measuring circuit at the time these points 
were taken, combined with the fact that the 
deflections were small because small tube currents 
had to be used at the high voltages. The differ. 
ence in the ordinates of the two curves is the 
specific intensity of the La doublet radiation from 
the thin target; it is plotted as curve I in Fig, 2. 


Ill. THE CORRECTIONS 
Fluorescence 


It is necessary to make allowance for the silver 
La radiation that arises from fluorescence of the 
target foil under the action of the x-ray continuum 
emanating from the beryllium backing. Let 
I(vo, v)dv be the power per unit solid angle 
emitted by the beryllium disk in the frequency 
range between v and v+dy, when ¥ is the high 
frequency limit of the continuum. Then the total 
fluorescent power emitted as La radiation by the 
silver foil is 


vo ee 2 aby Liit 
ure f avf ao f ——-I(v, v) 
vL 0 0 pv 


X {1—exp (—p~Xo sec 4)} sin 648, 


where uz4= the fluorescence yield of silver for its 
Lae radiation; 4=the absorption coefficient of 
silver ; uLz11 = the part of the absorption coefficient 
of silver due to absorption by the Ly shell; 
viy11 = the frequency of the silver Ly absorption 
edge ; Xo=the thickness of the silver foil ; @=the 
angle between an x-ray and the normal to the 
silver foil; ¢=the azimuth measured about the 
normal to the foil. 

If the crystal subtends at the target a solid 
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angle AQ, then the fluorescent power striking the 
crystal is 
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vIn MY 








«/2 
f {1—exp (—uXo sec 6)} sin 6dédv. 
0 


I(vo, v)AvAQ is the power striking the crystal in 
the range Av in the continuum. (A small and 
easily evaluated correction must be made for the 
absorption of the continuum in the silver.) This 
power has been measured, in arbitrary units, for 


the frequency corresponding to 4.3A. Let », be - 


this frequency, and J, be the power per unit solid 
angle per unit frequency interval at v, (corrected 
for absorption in the silver). The arbitrary units 
used to measure 7,Av and the uncorrected power 
of the La radiation are the same. Therefore the 
fluorescence radiation contributed to the ob- 
served specific intensity an amount 


vo MLUIVLIII 

$(r)= Sure f” “1(00,») 1 -GuXe) de, 
Y¥Lin “MY 

where 


«/2 
G(uX 0) -{ exp (—uXo sec 8) sin 6d8@. 
0 


The function G(y) has been tabulated by Gold." 

The data on fluorescence yields cited by 
Compton and Allison,” in conjunction with the 
relative intensity measurements of Parratt,? indi- 
cate that u#z¢=0.17. Van Dyke and Lindsay™ 
have determined vz for silver, and the work of 
Andrews" makes available the values of u and 
uly in the pertinent spectral region. The spectral 
band width admitted to the ionization chamber at 
ve was Av=3.7 X10" sec... 

Kulenkampff's'’ conclusion that I(vo, v) is pro- 
portional to v»— v is based on data corrected for 
absorption in the target. It implies that the 
intensity at a given frequency in the continuum 


1 E. Gold, Proc. Roy. Soc. A82, 62 (1908). 
* A. H. Compton and S. K. Allison, X-rays in Theory and 
_— (D. Van Nostrand, New York, 1935), pp. 488, 


e S. D. Van Dyke and G. A. Lindsay, Phys. Rev. 30, 562 
“C.L. Andrews, Thesis, Cornell, 1938. 
“© H. Kulenkampff, Ann. d. Physik 69, 548 (1922). 
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varies linearly with tube voltage. Measurements 
in the present investigation show that the in- 
tensity at »., measured at an angle of 85° with 
the cathode rays, increases less rapidly than the 
tube voltage up to 20 kv and is constant as the 
voltage ranges from 20 to 30 kv. Kulenkampff's 
formula for I(vo, v) is therefore not applicable to 
this case because of absorption of the continuum 
in the beryllium. Moreover, Kulenkampff worked 
only with the radiation emitted in a direction 
perpendicular to the cathode rays, while the 
fluorescence is excited by radiation passing 
through the foil in all directions. The numerical 
calculations show, however, that the foil absorbs 
strongly only those rays passing through the 
target rather obliquely. For this reason, and be- 
cause the cathode rays in the beryllium experi- 
ence considerable diffusion, it seems safe to ignore 
the dependence of I(vo, v) on @; the best expres- 
sion for I(vo, v) to use on the basis of available in- 
formation is I(vo, v)=(vo—v)I./(vo—v-). The 
function ¢(vo) was evaluated by plotting the 
integrand of (1) for several values of vo and then 
integrating with a planimeter. The results appear 
in column 3 of Table I. The quantity U in the 
table is the ratio of the tube voltage to the Ly 
ionization potential, while ¢(vo) is expressed in 
cm of galvanometer deflection per ma of tube 
current. 


Diffusion 
The fact that slow cathode rays are likely to 
experience large deflections in traversing the 


silver film makes the effective film thickness 
uncertain at the lowest voltages. For potentials 


TABLE I. Observed L doublet intensities, corrections, 
reduced intensities, and cross sections. 











U Tobs (v0) Xo/X Cs cr io Q 
cm/ma cm/ma cm/ma 
1.5 5.5 0.1 0.58 1.39 0.98 4.3 0.46 
2 11.8 0.3 0.62 1.05 0.96 voll 0.77 
2.5 14.0 0.4 0.65 1.01 0.95 8.5 0.91 
3 14.7 0.6 0.70 1.00 0.95 9.4 1.00 
3.5 14.8 0.7 0.74 1,00 0.94 9.8 1.04 
4 14.2 0.8 0.77 1.00 0.93 9.6 1.02 
5 11.5 1.0 0.83 1.00 0.93 8.1 0.86 
6 9.9 1.2 0.87 1.00 0.93 7.0 0.75 
7 9.2 1.2 0.90 1.00 0.93 6.7 0.72 
8 8.8 1.2 0.92 1.00 0.93 6.5 0.69 
9 8.5 1.3 0.94 1.00 0.92 6.2 0.66 
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above about 10 kv (U=3.0) a satisfactory calcu- 
lation of the effective thickness of the foil can be 
made. Throughout the discussion of this and the 
succeeding calculations the notation of Webster, 
Hansen, and Duveneck* will be followed as 
closely as possible. 

According to Bothe,'* the most probable deflec- 
tion experienced by a cathode ray in traversing a 
thickness x of silver is \=18.5U-'x!, where X is 
in radians, and x is in microns. Let 0 be the 
deflection of any particular cathode ray at a 
depth x in the film. For this ray, the path length 
in a thickness dx at depth x will be dx sec 0. The 
effective thickness of the target will therefore be 


Xo 
X= (sec O)xdx. 


0 


Here X» is the thickness of the film and (sec ©), 
is the value of sec 0, at a particular depth x, 


averaged over all of the cathode rays. For the. 


distribution-in-angle of the deflected cathode 
rays, Bothe'® gives 


n 
nd = exp (— @2/2n2)dQ, 


Tv 


where ndQ is the number of particles scattered 
into an element of solid angle dQ when mp elec- 
trons impinge normally on the film. This distri- 
bution function leads to a physically absurd 
result unless © remains small. Neher’s'’ data 
indicate that the number of electrons emerging 
along paths nearly parallel to the foil surface is 
small, and that 


ndQ « exp (— ©@?/2A*) cos OdQ 
is more nearly correct. Then 


(sec O) ay 


w/2 
f sec 0 exp (— 0?/2A*) cos @ sin OdO 
0 





x/2 
f exp (— 07/2?) cos © sin OdO 
0 


16W. Bothe, Handbuch der Physik XXII/2 (Springer, 


Berlin, 1933). 
17H. V. Neher, Phys. Rev. 37, 655 (1931). 
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This reduces to (sec 9), =2B(p)/A(p), where 


r/2p 
B(p) -{ exp (—w?*) sin pudu, 


r/2p 
A(p) -{ exp (— 1?) sin 2pudu 


and p=Av2. 
For small values of p, B(p) and A(p) are ap- 
proximately equal to B,(p) and A.(p), where 


B.(p) -f exp (—w?) sin pudu 
p/2 
=exp (—2°/4) [exp (ss 


and A..(p) =B,(2p). The function fo”? exp (s?)ds 
has been tabulated by Dawson.'® Investigation 
shows that B,(p) and A.(p) approximate B(p) 
and A (p) within 0.1 percent as long as p=0.5, but 
that the approximations become rapidly worse 
for larger values of p. In the interval 0.6=p=2.0, 
the Gauss 4-point numerical method of inte- 
gration’® yields values of B(p) and A(p) that are 
in error by less than 1 percent. 
The effective thickness of the target is 


Xo 


X=24 [B(p)/A(p) lax 


0 
bk 4.22/U 


=— B A > 
170 Je PL B(p)/A (p) lap; 


X was evaluated for various values of U by 
further application of the Gauss 4-point method, 
after determining pB(p)/A(p) for the appropriate 
values of p. The results appear in the fourth 
column of Table I. 


Retardation 


Let i(U)dx be the La doublet intensity from an 
ideally thin target of thickness dx. A real target 
of thickness X will give an intensity (i(U’))», 
where U’e is the energy of the electron after it has 
penetrated to a depth x, and (¢), denotes an 
average of i taken from x=0 to x=X. To a first 


18H, G. Dawson, Proc. Math. Soc. London 29{IT], 519 


(1897/98). ; 
19 K, Mader, Handbuch der Physik, ITI (Springer, Berlin, 


1928). 
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approximation,® 


(U)44X di sdU" 
i(U) w= i(U) + —(—) 
(i(U"))m WON ax Jyiw 

The observed intensities, corrected for fluores- 
cence and diffusion of the cathode rays, can be 
represented toa sufficient degree of approximation 
by the analytical expression 1(U)=KU~ log U. 
Therefore (i(U’))~=2(U)+g(U), where 


dU’ 
o(U) =4KXU-*(1 —log u(—) 
dx Ula 


and the observed data, after correction for the 
effects of fluorescence and diffusion, can be 
further corrected to compensate for retardation by 
multiplying them by c,=i(U)/[1(U)+g(U)]. 
Terrill?° and Williams** have found that, for 
cathode rays of about the energy used in the 
present investigation, (dU’/dx) « (1/U’)™. Terrill 
finds m=1, while Williams gives m=0.9, and re- 
ports a constant of proportionality much smaller 
than Terrill’s. Since Terrill worked with silver 
and other metal foils, and Williams’ work was 
done in a cloud chamber, it seems better to use 
Terrill’s results here. Calculation on this basis 
yields the values of c, shown in Table I. 


Rediffusion 

When cathode rays fall normally on a solid 
body, some of them, experiencing within the 
body deflections greater than 90°, are ‘‘rediffused”’ 
so that they appear to have been reflected from 
the body. The presence of the beryllium backing 
therefore causes some of the cathode rays to pass 
through the target twice. The fraction p of 
electrons experiencing rediffusion when cathode 


* H. M. Terrill, Phys. Rev. 22, 101 (1923). 
*E. J. Williams, Proc. Roy. Soc. A130, 310 (1930-31). 
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rays are incident normally on a plane surface of a 
material is called the ‘“‘rediffusion constant” 
of the material. Neher’s'? measurements on 
beryllium show that for 130-kv_ electrons, 
p=0.0248, and for 70-kv electrons, p=0.0291. 
Using data obtained by Schonland” for heavier 
elements and voltages from 10 to 100 kv, 
Webster, Clark, and Hansen? had previously 
estimated that p=0.043. A reasonable estimate, 
therefore, is that in the range from 10 to 30 kv 
p=0.035 for beryllium. Following Webster, 
Clark, and Hansen, let i(U) =io(U)+r(U), where 
to( U) is the x-ray intensity per unit length of path 
of the cathode rays, and r(U) is the intensity per 
unit thickness of target film. (The diffusion of 
rediffused electrons is a second-order effect which 
will be ignored.) It follows? that 


rU)=29 f F(W)-i9( UW) dW, 
1/U 


pF(W) being the fraction of incident electrons 
which are rediffused with energies between UeW 
and Ue(W+dW). All of the available evidence? 
on the energy distribution of rediffused electrons 
indicates that F(W) = (¢+1)(¢+2)(1—W) W?isa 
suitable approximation if g is an appropriate 
integer. For beryllium g=3 seems to be a good 
choice. An estimate of r(U) satisfactory for the 
present purpose can be obtained by using for 
io UW) the empirical approximation io(U) 
=KU- log U. Setting ¢,=io(U)/[io(U)+r(U)], 
so that i9(U)=c,-i(U), one obtains for c, the 
values shown in Table I. 


IV. RESULTS 


Curve I in Fig. 2 shows the observed specific 
intensity of the La doublet (corrected for the 
continuous background). To the ordinates must 
be applied first the correction for fluorescence, 
(vo), and then the corrections for diffusion, 
rediffusion, and retardation. Curve II in Fig. 2 is 
a plot of the observed La intensity minus the 
intensity due to fluorescence. The ordinates of 
curve III, the intensities per unit effective target 
thickness after all corrections have been made, 
are obtained from those of curve II by multiplying 
them by c,-c,-Xo/X. The finally corrected re- 


= B. F. J. Schonland, Proc. Roy. Soc. Al08, 187 (1925). 
























sults, which are proportional to the ionization 
cross section Q, are listed under ip in Table I and 
are replotted in Fig. 3, the ordinates here being 
arbitrary units chosen to make Q=1 when U=3. 
The theoretical values of ionization cross section 
calculated by Burhop*‘ are shown as crosses ; they 
are fitted to the experimental curve at U=3. 
Burhop states that his results have been evaluated 
numerically to ‘‘a few percent.’”’ The measured 
values of Q are tabulated in Table I. 

There are four sources of error in the experi- 
mental results, namely : 


(1) the measurement of target potential, 
(2) the measurement of tube current, 

(3) the measurement of x-ray intensity, and 
(4) the corrections. 


The target potential was measured with a low 
voltage laboratory standard voltmeter and wire- 
wound multiplier. Checks with a potentiometer 
and Wheatstone bridge showed the calibration of 
the combination to be correct within 0.25 percent, 
when account is taken of the corona and surface 
leakage in the multiplier. Errors in the measure- 
ment of target potential therefore make no 
significant contribution to the experimental error. 

The target currents could be read within about 
1 percent. 

Theelectrometer tube circuit introduced uncer- 
tainty into the intensity measurements through 
fluctuations and through changes in sensitivity. 
Frequent checks on the sensitivity ensured that 
the errors due to changes in sensitivity are cer- 
tainly less than 2 percent. Fluctuations are par- 
tially neutralized by drawing a smooth curve 
among the points ; they render the location of the 
curve uncertain by not more than 2 percent. This 
estimate includes an allowance for zero drift. 

The geometrical thickness Xo of the target is 
uncertain by about 15 percent, but it transpires 
that a 15 percent error in Xo causes an error of 
little more than 1 percent in the measured value 
of ionization cross section at U=9 relative to 
that at U=3. The errors introduced at other 
values of U between 1.5 and 9 are even smaller. 
Departures from the assumed law of distribution- 

in-angle of the diffused cathode rays may intro- 
duce further errors, but no estimate of their size 
can be made except that they are probably not as 
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large as 1 percent. They will therefore be 
neglected. 

The corrections for rediffusion are so small and 
are based on such accurate information that they 


do not contribute significantly to the experj. © 


mental error. The same is true of the corrections 
for retardation, except for that at U=1.5, which 
renders Q at that voltage uncertain by about 
5 percent. 

The square root of the sum of the squares of the 
errors that have been mentioned is 6 percent at 
U=1.5 and 3 percent for the larger values of U. 
These are the experimental uncertainties in the 
ordinates of Fig. 3 at small values of U; at large 
values of U there is in addition the possibility of 
errors introduced by the assumptions underlying 
the correction for fluorescence. Concerning these 
errors, all that can be said is that they are 
negligible for small values of U, say U<7, and 
that the computed correction for fluorescence is, 
if anything, too small. (It is perhaps worth 
mentioning here that Burhop’s results were not 
consulted until after the calculation of the cor- 
rections had been completed.) For U>7, the 
curve does not fall off as rapidly as the theory 
predicts ; a 50 percent increase in the correction 
for fluorescence would be necessary in order to 
bring the experimental point at U=9 into 
agreement with the theoretical prediction. With- 
out making actual measurements on the intensity 
of the x-ray continuum from beryllium in direc- 
tions making various angles with the target, over 
the range of tube voltages used in the present 
experiment, one cannot decide with assurance 
that, for U=9, the present correction for fluores- 
cence is in error by less than 50 percent. It should 
be pointed out, however, that the experimental 
results of Webster, Hansen, and Duveneck,? of 
Smith,” of Tate and Smith,™* of Webster, 
Pockman, and Kirkpatrick,?® and of Webster, 
Pockman, Harworth, and Kirkpatrick** all show 
that the theoretical calculations (based on the 
Born approximation) predict cross sections that 
decrease too rapidly with increasing U when-Uis 
large. 

In conclusion, it seems of interest to compare 


%P. T. Smith, Phys. Rev. 36, 1293 (1930). 

“J.T. Tate and P. T. Smith, Phys. Rev. 39, 270 (1932). 

2% Webster, Pockman, and Kirkpatrick, Phys. Rev. 4, 
130 (1933). 

26 Cited by Burhop, reference 4. 
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the present measurements on the Lin state of 
silver with those of Webster, Hansen, and 
Duveneck® on the K state of the same element. 
Two facts are clearly shown : 

(1) The cross sections for excitation to the K 
state and the Lym state are both greatest for 
cathode rays whose kinetic energies are about 3.5 
times the ionization energy of the state in 


question. ene 
(2) When U>3.5, the decrease in ionization 
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cross section with increasing cathode-ray energy 
is more rapid for the Lyn state than for the K 
state. 
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The invariant field theory of Part II is interpreted, in 
agreement with F. Bopp, as Maxwell's theory with a linear 
differential relation between the fields ZE, B and D, H 
involving a new constant k which measures the reciprocal 
radius of the electron. The former “‘mesonic field’ of 
minimum frequency vo=c/2 represents polarization of 
the vacuum. The electron is a singularity in the D, H 
field whereas FE, B remain finite. Instead of obeying 
dynamical equations of motion, the electron moves under 
the condition that the Lorentz force vanishes identically 
on the singularity, so that no work is done on the particle. 
All energy is located in the field. In this respect the theory 
is unitary. Electromagnetic and inert mass are identical. 


1. CLASSICAL FIELD THEORY 


HE modification of electrostatics proposed 

in Part I and its electromagnetic continu- 
ation discussed in Part I]! rest on the assump- 
tion that vacuum is polarizable, as described by 
linear differential relations between the vectors 
E, B and D, H (for details see Section 6): 


D=E-k"OE, H=B-k“OB, (1) 


where (J is the Laplace operator in x, y, 2, ict. 
The constant k of dimension [/-'] determines 
the minimum frequency »o=kc/2m of waves of 
polarization or ‘‘meson waves.”’ k also plays the 
role of the reciprocal electronic radius, although 
the charges ¢ are condensed in mathematical 


1A. Landé, Phys. Rev. 60, 121 (1941). A. Landé and 
L. H. Thomas 60, 541 (1941). 


In contrast to Dirac’s classical electron which is subject 
to advanced and retarded potentials and displays self- 
acceleration, the field theory works with retarded poten- 
tials only, and self-acceleration is avoided. Stable equi- 
librium between electrons and radiation is granted by 
spontaneous and induced transitions, similar to Einstein's 
derivation of Planck’s radiation formula. In spite of dis- 
playing a magnetic moment the electron does not have 
magnetic self-energy, so that its radius is the ordinary 
electrostatic radius 1/k=e/2mc*. In contrast to Born- 
Infeld’s non-linear theory, our field equations allow a 
Fourier representation as a basis for the quantum theory 
of Part IV. 


points only. The simplicity and naturalness of 
our approach are demonstrated by the fact that 
the same modification of Maxwell's theory has 
been proposed independently and simultaneously 
by F. Bopp.? Whereas we began in Part | with 
a Fourier representation of the field of a point 
particle with finite self-energy, Bopp started 
from a formal generalization of the Lagrangian 
function of the field E, B= f.g, namely, 


L=—(1/16) { (fas)? +k-*(9 fap/O%4)?}+Jaga (2) 


where J is the 4 current and ¢ is the 4 potential. 

The relation to other field theories (Maxwell, 
Born-Infeld) become obvious if the vectors E, H 
of Part II are called EZ, B, and the vectors E”’, 
H" are called D, H. Our “‘meson field” E’ = D—H 


2 F. Bopp, Ann. d. Physik 38, 345 (1940). 









176 





TABLE I, Notation. 


Our former E,H;V,A E”",H";V",A" E',H’; V',A’ 
Bopp E,B;¢,a@ D,H;¢,A F,G;V/k,U/k 


and H’=H-—B then represents the electric and 
magnetic polarization of the vacuum (see Table I). 
In contrast to Born-Infeld, our relation between 
D, Hand E, B is linear so that it is easy to obtain 
a Fourier representation (not discussed by Bopp) 
with subsequent quantization. Infinities are 
avoided automatically without resorting to arbi- 
trary cutting-off procedures. 

The new field theory is unitary insofar as all 
energy is supposed to be located in the field, with- 
out additional ‘‘mechanical”’ energy of particles. 
Instead of introducing mechanical equations of 
motion for the electrons, we postulate that the 
field E and B (which is responsible for Lorentz 
force and work on the particles) shall vanish at 
all times on the point singularities themselves. 
As a consequence of this condition we were able in 
Part II to derive a quasimechanical equation of 
motion: 


m(#)9= external force + (2€/3c*)(d%/dt)) (3) 


valid at a particular time (¢=0). m is an abbre- 
viation for ke?/2c? and represents the electromag- 
netic mass of the surrounding field. 

The simplest solution of the field equations is 
an electrostatic field of spherical symmetry 
surrounding a point charge e (II, Section 4): 


D=eé/r, E=e/r?—e(1/r?+k/r) exp (—hr). 


In contrast to D, the field E remains finite at 
r=0, namely, E=ek?/2. However, similar to 
Born-Infeld, the field component E, jumps from 
the finite value +|£| to —|E| when passing 
through the singularity, with the average value 
E,=0 on the singularity itself. This discon- 
tinuity preserves the individuality of the par- 
ticle. The electric field energy of the particle at 
rest has value W=ké/2=mc’. 

Two point charges €, and ¢2 at distance r have 
field energy 


W=khe?/2+ke?/2+(ee/r)-[1—exp (—kr) ]. (4) 
Two opposite charges +¢ and —e therefore yield 


W=mc?+mc— é/r for large r, (5) 
W=mckr for small r. 


Hence, when r decreases from « to 0, field 
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energy of value 2mc? is released. It thus turns 
out that transformation of mutual electrostatic 
into radiation energy with neutralization of +¢ 
and —e (annihilation) is possible already in the 
classical domain. 

Another significant feature of our theory 
concerns electrons vibrating with high fre. 
quencies v>vo=kc/2x. The energy emitted jn 
this case is reduced by a factor (II, Section 6) 


1 — (1 — 0?/v?)! (6) 


as compared with the normal Maxwell-Lorentz 
radiation. For v>vo' the reduction factor ap- 
proaches v?/2v?, so that energy is emitted at the 
rate vy? rather than yv* per second, or »v rather 
than v* per period. This is of importance for the 
quantum theory of stationary states which 
always presumes that the energy emission during 
a single period is small compared with the energy 
of the state itself. 

Our theory in its classical form does not dis- 
tinguish between self-field and external field, 
except for the special case of uniform motion, 
On the other hand, the quantum theory of radi- 
ation is dualistic. Here one considers electrons 
in uniform motion with definite self-energies 
W = W,(1—8?)-? surrounded by a “pure field” 
with constant Fourier amplitudes. Transitions 
between stationary states result from mutual 
perturbations between pure field and particles. 
The difficulties of the quantum theory of radia- 
tion arise primarily from the task of describing 
the interaction between external field and 
mechanical particles in a dualistic fashion, 
although the classical theory of the electron in 
the field is unitary. The difficulty is solved in 
Part IV by considering as zero approximation 
particles without rest mass which therefore 
travel with the velocity of light. When the inter- 
action with the field is introduced, stationary 
states of the particles with any velocity smaller 
than c become possible. These states may be 
obtained from states with zero velocity by 
Lorentz transformations, and the rest mass now 
has a definite value. Our special type of semi- 
unitary field theory in which the particles retain 
their individuality as singularities although all 
energy is located in the field, offers a satisfactory 
background from the classical and quantum 
points of view. 
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2. THE PROBLEM OF SELF-ACCELERATION 
For vanishing external field Eq. (3) reads 
dié/d(ct)=aé% with a=3mc?/2é. (7) 


The solution of this non-relativistic equation is 


dx/d (ct) =exp (act)+ 8 (7’) 
where 6 is the velocity ratio v/c att=—o. A 
free electron thus should be able to accelerate 
itself from any initial velocity v at :=— © to 


larger and larger velocities, the beginning phase 
of the process being described by (7’). Rela- 
tivistic corrections are needed for higher veloc- 
ities only. Dirac* has emphasized the funda- 
mental importance of self-acceleration in his 
relativistic theory of the classical point electron. 
There he assumes that the force on the point 
charge is determined by the difference of retarded 
and advanced potentials, in order to get rid of 
the infinite self-energy. As a consequence he 
obtains self-acceleration. However, the whole 
idea of a free particle acquiring larger and larger 
velocities on its own accord is so non-physical 
that any theory yielding self-acceleration might 
be rejected almost a priori. Yet, self-acceleration 
is an inevitable consequence of any theory in 
which the electron is subject to a differential 
equation of motion of finite order because the 
equation of motion at the beginning of the process 
for small velocities is always approximated by 
(7) with the solution (7’). It is a decisive ad- 
vantage of the present field theory that the 
singularity does not obey a differential equation 
of motion. The force on the electron at ¢ depends 
on the whole path before ¢ rather than on the 
motion in the immediate neighborhood of ¢. 

In II, Section 4 we discussed an electron 
moving with coordinate 


x=ut+$ft+éige (7”) 


during a short time ¢ near ¢=0 and with small 
coefficients u, f, and g whose meaning is u = (Z)o, 
f=(®)o, g=(d#/dt)y at t=0. The self-force was 
found to be 
self-force = — (e&k/2c*) - f+(2e/3c*) -g 
= —m(#)o+ (2€/3c*) (dz/dt)o 


where m is an abbreviation for (@k/2c*) repre- 


*P. A. M. Dirac, Proc. Roy. Soc. 167, 148 (1938). 
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senting the finite mass of the point charge «. 
Our requirement of a vanishing total force means 
vanishing self-force for a free electron, so that 
the latter at ‘=0 satisfies the equation 


m(#)o= (2€/3c) (d#/dt) . 


This equation does not permit conclusions 
about self-acceleration before or after ¢=0. 
Whether self-acceleration is possible depends on 
the question whether the self-field vanishes during 
an accelerated motion. Before answering this 
question of dynamics we must first find the cor- 
rect kinematic description of a relativistically 
accelerated motion. 

Let us assume tentatively that the motion of 
a certain free particle were actually of the form 
(7’) at least for large negative ¢ and for small 
values of the initial velocity 8. Let us ask also: 
What is the correct relativistic continuation of 
this accelerated motion? The answer may be 
found by relativistic transformations without 
resorting to dynamics. Since the electron is free, 
the acceleration at any time must depend only 
on the proper time s measured on the particle 
itself. For large negative ¢, however, ¢ and s are 
identical provided that B<1. Instead of (7’) we 
therefore may write relativistically 


dx/d(cs) =exp (acs)+8, 


d(ct)/d(cs) =1 for s=— @. (8) 


With reference to a Lorentz system xt in which 
the electron is at rest at the beginning of the 
process we also have 


dxo/d(cs) =exp (acs), 


(8’) 
d(cto)/d(cs) =1 for s=— &. 
Generalizing these equations tentatively for all 
s values we write 


dx/d(cs) = S(e**+6), d(ct)/d(cs)=C(e**+b), 
dxo/d(cs) = S(e***), d(cto)/d(cs) = C(e**), 


where 6(8) is a function of 8 which coincides with 
8 for small 8, whereas S(z) is a function of z which 
coincides with z for small z like a sine function, 
and C(z)—1 for small z like a cosine function. 
We now write down relativistic transformation 
formulae from xt to xoto 


dx = [dxo+d (cto) “Bi — p*)-, 
d(ct) = [dxo-B+d(cto) }(1 — B*)-1, 
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Dividing by d(cs) and writing z for e** we obtain 


S(s+b) =[S(z)+C(z) -8]-(1—6*)-4, 
C(2+b) =[S(2)8+C(z) ]-(1—6?)-4. 


These equations represent rotations in Minkow- 
ski space. Indeed, if we write 


(1—f*)-!=cos (ib) = Cosh (8), 


B(1—6')-+= —i sin (ib) =Sinh (b), Pann 0=8. 


we have 


S(z+b) = S(z) Cosh (b)+C(z) Sinh (8), 
C(z+5) = S(z) Sinh (6) +C(z) Cosh (8). 


These formulae are valid only if the functions S 
and C are Sinh and Cosh. The correct relativistic 
continuation of (7’) in the xt system, therefore, 
is 


dx /d(cs) =Sinh (e**+5), 
d(ct)/d(cs) =Cosh (e**+5), 


where B=v/c describes the original constant 
velocity of the electron in the xt system at 
t=— , and s is the proper time. Equation (9) 
represents the only accelerated motion possible 
for a free electron with (7’) as non-relativistic 
approximation for !~ — «. As a consequence we 
obtain the following invariant description of 
self-acceleration : 


d*x,/d(cs)?=a-exp (acs), (9’) 


where x, is the Lorentz system in which the 
electron is momentarily at rest. 

Next let us calculate the self-force according 
to the field theory if the electron is relativistically 
accelerated according to (9). In order to simplify 
the mathematical problem let us discuss only the 
initial part of the process between ¢= — © and 
t=(1/ac)-In C where C is a small positive 
number. During this initial phase we may use 
the approximation (7’). Shifting the zero point 
of the time scale we may consider the motion 


dt/d(cr) =C exp (acr), (9””) 
§=(C/a)Lexp (acr) —1] 


during — © <7<0. The zero point of the & axis 
coincides with the position of the electron at 
7=0. Furthermore, instead of discussing the 
self-force during the whole time 7 <0, we cal- 
culate it at r=0 itself. In the original time scale 
this is the time ¢=(1/ac) In C where C is any 


b=Tanh-' B, (9) 
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small number, so that what applies to !=0 also 
applies to a continuity of times ¢ in the original 
time scale. 

In order to find the self-force of the motion 
(9) at r=0 we apply the method of retarded 
potentials (Part II, Section 4) and proceed jp 
three steps. 

(1) If the exponential function in (9’) jg 
expanded into a series we have 


&=(C/a)[acr+}(acr)*+§(acr)?]. 


If we omit all higher terms, the result was 
already found in II, Section 4: 


Self-force = — (@k/2c*) (E) 0+ (2e/3c*) (dé/dr), 
= — &a®C[ (k/2a) — (2/3) ]. 


The self-force vanishes for a =0 (uniform motion) 
and also for a=3k/4 (self-accelerated motion), 
that is, for a=3mc?/2é which is Dirac’s value. 

(2) If the next expansion term (acr)*/24 js 
added to the series, the self-force becomes infinite 
for every value of a except a=0. 

(3) If the complete exponential function 
exp (acr) is used the self-force becomes finite 
again due to the fact that successive expansion 
terms contribute infinite terms of opposite signs, 
The resulting self-force according to the method 
II, Section 4 is 


—eéC- f “7 (RR) { (aR“e-**) + (a Re) 
— (Ra) + (R“e-*®) | d(kR) 

or after integration: 

Self-force = 

(a?+k?)'—a (a?+k*)i—k* 1 


a 3a’ 3 


—éa®C 








If we had used the exact relativistic motion (9) 
rather than (9’), the bracket would contain an 
additional term linear in the small quantity C. 
For small-but finite C the self-force vanishes only 
if the bracket is zero. The latter, however, does 
not vanish for a=3k/4 nor for any other value 
of a save a=0 (or a<0 if the small correction 
term linear in C is added). That is, our field 
theory does not allow self-acceleration in its 
beginning phase. Therefore, the whole process 
of self-acceleration is not compatible with the 
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field theory. It will be noticed that the bracket 
in the last equation for the self-force would 
vanish for a=3k/4, if the positive root (a?+k?)! 
were replaced by the negative root. However, 
because of our application of retarded potentials 
the root appears with positive sign, and we arrive 
at the gratifying result that self-acceleration is 
impossible in our field theory. 


3. NEGATIVE ENERGY DENSITY 


It seems impossible to devise a field theory of 
charged particles without violating the require- 
ment that the density of the field energy ought 
to be larger in the presence of a field than with- 
out field, that is, be positive definite in the 
presence of a field. As an example we mention 
Born-Infeld’s field theory which leads to negative 
energy densities‘ wherever the vectors D and B 
have certain large values. (In our theory negative 
field energy is connected with large frequencies 
y>vo=kc/2x rather than with large field inten- 
sities). Dirac* in his classical theory of the elec- 
tron simply assumes that an infinite amount of 
negative energy of unknown origin is condensed 
near the point electron so as to counterbalance 
somehow the infinite positive Maxwell energy. 
He thereby denies that the mass of the electron is 
of electromagnetic origin. 

The negative energy density of the meson 
field in our theory seems to lead to certain dif- 
ficulties. Consider an electron vibrating with a 
frequency »v larger than vo for a limited time only. 
As shown in Part II, it then emits a group of 
Maxwell waves E,B and a group of meson 
waves of polarization D—E, H—B. The two 
groups travel with different group velocities and 
become separated at some distance from the 
source. The Maxwell group carries more energy 
than was emitted by the electron; the meson 
group with its negative energy re-establishes the 
energy balance. A resonator exposed to the 
Maxwell group alone then might transform the 
Maxwell energy into mechanical energy of the 
resonator, thereby reducing the energy remaining 
in the field to larger and larger negative values. 
This objection does not hold, however, in our 


*M. Born, Proc. Roy. Soc. 143, 410 (1934). See also 
W. Heitler, Quantum Theory of Radiation (Oxford Univer- 
sity Press, 1936), Eq. (5), p. 237. 
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unitary theory where all energy is located and 
conserved in the field, and the work on the 
resonating charged particle is zero. From the 
energetic point of view our position with respect 
to negative field energy is the same as that of an 
observer above sea level who discovers that 
weights may be dropped below zero level at the 
expense of other weights raised above sea level. 
Difficulties appear, however, if one tries to 
derive a statistical distribution of photons and 
mesons over positive and negative levels, at 
least if he wants to apply the same statistical 
methods which had been developed for positive 
levels. 

In case of Boltzmann statistics the relative 
number of particles on two energy levels is 


fi1/fe=exp (W2- W,/kT), 


and remains unchanged if both energies are 
counted from a new zero point. Bose statistics, 
however, yield the following formula for the 
average number of particles of energy W: 


n=1/[exp (W/kT)—1), 


where W is an absolute energy value without an 
additional constant. 7;/f/2 depends on the ab- 
solute position of the zero energy. Furthermore, 
Bose’s formula is meaningless for negative values 
of W unless a new interpretation ofa “negative 
number of particles’’ is invented. However, one 
must remember that Bose particles are “par- 
ticles’ only in a very restricted sense, and that 
the success of Bose’s procedure is more or ‘ess 
accidental. Bose showed that the same Planck 
radiation law which was derived from an equi- 
librium between radiation and radiating matter 
(Planck’s and Einstein’s derivations) could also 
be obtained from certain statistics applied to 
“particles without individuality.’’ We cannot ex- 
pect that a modified Planck law (asked for by 
modified field theories) should still be derivable 
from Bose’s corpuscular method. This applies in 
particular to our case where ‘‘pure field’? and 
“electronic field’’ are undistinguishable within 
the range 1/k of the electronic “‘radius.’’ Radia- 
tion of wave-length \=22/k may be treated in 
connection with electrons responsible for the 
radiation equilibrium, similar to Einstein's 
derivation of Planck’s law, as follows. 

















4. MODIFIED EINSTEIN EQUILIBRIUM 


An electron vibrating with frequency v> vo 
=kc/2m emits negative meson energy at the 
ratio of (1—»?/»*)! to 1 in comparison to the 
normal Maxwell radiation, according to (6). 
This ratio has another significant meaning. The 
number of Jeans proper vibrations per unit cube 


and per wave-length interval dd is 
dZ = (8x/X*)d(1/d). 
Since c/A= v for Maxwell waves, and 
c/A= (v®— v7)! 


for meson waves, we obtain the following num- 
bers of levels within the same dy: 


dZ" =82rv'dv/c* for photons hy (10) 
dZ' = (8rv*dv/c*) (1 — v9?/v?)* for meaons (—hyv). 


dZ’ and dZ” turn out to have the same ratio as 
the numbers of mesons and photons emitted 
“spontaneously”’ by the vibrating electron during 
a certain time. Hence, if a group of electrons has 
emitted just one photon onto every one of the 
dZ”" levels within dv, the same group of electrons 
has also emitted just one meson onto every one 
of the dZ’ meson levels within dv. The total 
number of photons and mesons emitted is pro- 
portional to (dZ’’+dZ’) whereas the total energy 
emitted is proportional to hv-(dZ’’—dZ’). 


After this preparation let us discuss Einstein’s, 


derivation of Planck’s radiation law, modified by 
our field theory. Two electronic energy levels 
W, and W.<W, may be occupied by N, and 
N, electrons, respectively. Electronic transitions 
between the two levels produce emissions and 
absorptions of quanta +hv and —hy toand from 
the Jeans radiation levels within a certain re- 
sonance interval dv containing dZ” photon 
levels and dZ’ meson levels. An electronic transi- 
tion from W, to the higher level W, corresponds 
to an increase of the electronic amplitude of 
vibration and to a decrease of the Maxwell field 
strength | £’’| (=absorption of photons) and to 
an increase of the meson field strength | £’| 
(=emission of mesons). In terms of statistics 
this means that the probability of an electronic 
transition from W, to W, ought to be propor- 
tional to the difference of the numbers n’’dZ”’ of 
photons and the number n’dZ’ of mesons ready 
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to be absorbed from the resonance interval dy: 


Pra= N,[n"'dZ" nine n'dZ’ |. 


Opposite transitions a—d induced by the radia. 
tion will be proportional to the same difference, 
whereas spontaneous transitions will be propor. 
tional to 1-dZ’’—1-dZ’ [see (6) and (10)] go 
that 


Pus = Na[(n"’+1)dZ” —(n'+1)dZ']. 


In case of equilibrium we must have P=P,, 
in the average, that is, 


(a +1)dZ" — (ni +1)dZ" 


N./Na= 
n'dZ" —n'dZ' 





The bars indicate average numbers of photons 
and mesons per level. If equilibrium shall apply 
to photon processes alone and to meson processes 
alone we must have 


No/ Na= (a +1)/0" = (n' +1)/7’, 


hence, 7’’=7'’. That is, the average number of 
photons and mesons per Jeans level must be 
equal. If the electrons are subject to Boltzmann 
statistics N=const. exp (— E/kT) we obtain the 
following value for the average number of 
photons and mesons per level: 


n”’ =n’ =[exp (hv/kT)—1}". (10’) 
The total radiation energy within dv has density 
wav =n" hvdZ" +n'(—hv)dZ' = (w," —w,')dv 
8rhv*dv 1 
- Cc exp (hv/kT) —1 
X[1—(1—v0?/r*)*].  (10") 





For v>vo=kc/2e Planck's radiation formula is 
multiplied by the same factor (6) which already 
occurs in the classical energy emission of an 
electron for v> vo. The result is obtained without 
applying statistics to the pure radiation. 


5. STABILITY OF EQUILIBRIUM 


In order to show that the equilibrium is stable 
consider small deviations from the average 
numbers 
n'=n'+8', 

Nz == N, —A 


n"’ = n’ +6", 


N.=N.+4, 
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which are related by the condition 
A=Z's' —Z"5" 


where Z’ and Z” are the number of radiation 
levels involved (formerly called dZ’ and dZ’’). 
Photonic processes lead to an increase of 
the number of photons, d(n’’Z’’)/dt=P.a— Pra 
whereas dZ’ =0, that is, 


d(n"Z")/dt=(Nat+A)(0" +8" 412” 
—(N,—A) (a +8")Z". 


The terms of zero order in the deviations cancel, 
and the second-order terms shall be neglected. 
This leaves the first-order terms: 


d(8"Z")/dt= —(8"2Z")(N,— Na) +4- (20 +1)2”. 


Mesonic processes yield similarly d(n’Z’)/dt 
=P,,—P.» and dZ” =0, that is, after removing 
second- and zero-order terms: 


d(8’'Z’)/dt= —(8'Z')(Ne—N.)+A(2n'+1)Z’. 


At last, the electronic number JN, increases with 
probability Pi.—P.a, that is, after removing 
second- and zero-order terms: 


dA/dt= —A\(2n%" +1)Z”" —(2n’+1)2'} 
+(N,—N.)(6"2Z" —8'2’). 


The last bracket is —A. Writing C for the 
positive constant N,—N., we thus obtain 


(a) dA/dt= —A{C+(2a"+1)Z”" —(2%'+1)Z'} 
together with the former results 


(8) ds” /dt= —8"C+A-(2n"+1), 
(y) ds’ /dt= —8'C+A-(2n' +1). 


The factor of —A in (a) is positive since C>0, 
and ’’’ =n’, as well as Z’’ >Z’. Thus the absolute 
value of A will always decrease. The terms 
with A in (8) and (y) will therefore shrink to 
zero after some time, and from thereon the ab- 
solute values of 6’ and 6’ will always decrease. 
A small deviation from the average values 7’’, 
n’, N., Ny thus tends to disappear automatically 
according to the transition probabilities accepted 
before, and the equilibrium turns out to be stable. 


6. FIELD EQUATIONS 


We now return to our former notation E, H 
for Maxwell’s field E, B, and E”, H” for Max- 
well’s D,H (Table I). Both fields shall be 
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derived from potentials V,A and V”,A”’, re- 
spectively : 
E=-VV-A/c, 


H=VXA; (11) 
E" = —VV"—A"' Je; 
H"”=VxXA". (11 ) 


However, only the potential V’, A” shall always 
obey the Lorentz condition: 


(V-A)+V/c=R, (12) 
(V-A")+V"/c=0 (12’) 


where R on the right of (12) may be any scalar 
function of xyzt, to be restricted later. In our 
applications we usually consider the special case 
R=0 only. As a consequence of (11), (12) one 
obtains field equations for E,H and E”,H": 


VXE+H/c=0, (13) 
(V-H)=0; 
VXE"+H"/c=0, (13’) 
(V-H’’) =0; 
VXH-E/c=—-DOA+0R, (14) 
(V-E)=—-OV-R/c; 
VXH" —E"/c=—-OA", (14’) 


(V-E”) = — oO yr 


The right-hand sides of (14’) represent 4 times 
the current j/c and density p of the true charge 
(which later on is supposed to be condensed to 
world lines only). The right-hand sides of (14) 
represent 4m times the free current and density. 
OA and OV have the dimensions of potentials 
divided by the square of a length. They may be 
written in the form 


—-DA=kA’, (15) 
—-OV=kRV’; 

— DA" =4xj/c, (15’) 
=f} V" =4np; 


thereby defining a new “potential’’ A’, V’. 
Hence (k?A’+VR) is 4x times the free current 
density, and (k?V’—R/c) is 4” the free charge 
density. As a consequence of (14) (14’) we obtain 
continuity equations for free and true charge: 


(V-A’)+V'/c=—OR/F, (16) 
(V-7)+6=0. (16’) 
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From the potentials A’, V’ we may also derive 
a new “field” E’, H’: 


E'=-VV'—A'/c, H'’=VXA'. (17) 


(13) and (14’) are Maxwell’s equations. It has 
been assumed that both fields E, H(=E, B) 
and E”, H’’(=D,H) are derived from poten- 
tials, and that V’, A” satisfy the Lorentz con- 
dition. 

We now postulate a new relation between the 
two potentials, namely, 


V"=V-OV/R, A” =A-—OA/R* (18) 
from which follows 
E"=E-ODOE/P, H”=H-OH/k* (19) 


or D=E-—OE/k? and H=B-—OB/k* in the 
usual notation. Eliminating 0A and OV from 
(15) and (18) we obtain 


A=A"—A’, V=V"-V’. (20) 


Hence (12) must be the difference of (12’) and 
(16). That is, the right-hand side R of (12) must 
satisfy the condition 


—OR+FR=0. (21) 
From (20) we learn that 
E=E"—E’, H=H"—-—H’' (22) 


or E=D—-E’ and B=H-—H’ in the usual 
notation, which shows that the fields E’ and H’ 
determine an electric (P) and magnetic (M) 
polarization of the vacuum 


E'=4nrP, H'=—4nrM. (22’) 
Equation (16) now reduces to 
(V-A’)+V'/c=—R. (23) 


The field equations for E’ and H’ according to 
(17) and (23) read 


VX E'+H'/c=0, (24) 
(V-H’)=0; 


Vv XH’ — E'/c=—OA’—-VR, (24’) 
(V-E))=-OVW'+R/c. 


Subtraction of (15) from (15’) yields, because of 
(20), 


—DA'+k?A’ =4xj/c, (25) 
—OV’+k?V’ =4rp. 
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Equation (25) together with (15’) 


—OA"=4nj/c, —OV"=4rp _  (25’) 


shows that the true charge is the common 
source of the two otherwise independent fields 
E”,H" and E’, H’. True charge shall be con. 
densed on singular world lines only. 

The absolute values of the potentials A’, y’ 
have a physical meaning without additional 
constants since they occur in the field equations 
(14) whose right-hand sides read 


RA'+VR and #V’—R/c. (26) 


They represent 4x times the free current and 
free charge density. 

Since the fields E=E"’—E’ and H=H" —}’ 
determine the Lorentz force and work, the energy- 
momentum tensor ZJ must be the difference 
T=T" —T". In particular, the density of energy 
and flux are 


S=S"-S’, (27) 


w=w'— w’, 
where 


w’ = (1/4) {}(E°+H?+hA?+kV") 
+RV'/c—RV'/c+3R*}, (28) 


S' = (c/4) { LE’ XH']+kV’A'+ V'VR+A’'R/c}, 


and 
w’’ =(1/49r)}(E’'"+H"”), (28’) 
S” =(c/4r)[E” XH"). 


Equation (28) contains the right-hand side R of 
(12) and is more general than the expressions 
for w’ and S’ of Part II where we considered 
the case of R=0 only. 

Meson waves in vacuum may have longitudinal 
field components. Indeed, consider the special 
case 

A,’ =a,’ sin (2rvt+2xx/d), A,’ =0, 
A, =0, V’=v'-sin (2evt+2xx/)X). 


The amplitudes a,’ and V’ are not restricted by 
(12). The field is 


E,! = —(az'v/c+v'/d)2x cos (24vt+2rx/d), 
E,' =E,' =H,' =H,' =H,’ =0, 


representing a longitudinal electric wave. In 
case of Maxwell waves we would have a,” = —v", 
and E,’’=0. Maxwell waves do not have longi- 
tudinal components for two reasons: first, 
because A\=c/v, and second, because of the 
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Lorentz condition (12’). Meson waves do not 
have to satisfy the Lorentz condition, and their 
relation between ) and » is 


(v/c)?=1/d?+ (ke/2x)?. (29) 


7, ELECTRIC POLE AND MAGNETIC DIPOLE 


In Part II we discussed the potential and field 
of an electric point charge ¢. The field equations 
allow the following solution: 


Vv=e/r, V'=(e/r) exp (—kr), (30) 
E,"=¢/r, E,’=e(1/r+k/r) exp (—hr). 


At large distance this is the ordinary Coulomb 
field of a point charge e. In a similar way we 
now discuss the field of a magnetic dipole of 
moment yu. The field equations for R=0 allow 


the solution: 
dsl 

soon acer) seer) 
oz\r 


° ] 
A,’ =0, A,’ =y—(e*"/r), 
Oz 


Fi 
A,=—yp—(e*"/r) (31) 
oy 


with the magnetic fields 
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For large r this is the ordinary Coulomb field of 
a dipole of moment yu parallel x. Although the 
field H=H’—H’ does not vanish for r=0, the 
Lorentz force w@H/dx vanishes for r=0. The 





classical field theory remains unitary if we 
postulate as before that the electron moves so 
that the total Lorentz force and work on the 
singularity are always zero. 

In order to find the energy of the magnetic field 
we use (27), (28) for the energy density w=w”’ 
—w’ with function R=0. Integration over space 
yields 

W mag = w*h*}. 


The field energy of the point dipole is finite due 
to magnetic polarization of the vacuum. 

If the electron really had a magnetic self- 
energy like this, the common electric and mag- 
netic radius 1/k at which Coulomb’s law breaks 
down would be determined by the equation 


me=tek+ $urk'. 

Substituting the quantum values 
p=eh/2mc, @=ahc, a=1/137, (32) 
this would yield the following equation for 


5=2mc?/ ek: 


6 
—=}+4(ad)~* 
2 







solved by 6= 16.47; that is, 1/k = 16.47 (€/2mc). 
The radius would be more than sixteen times the 
electric radius, and the magnetic energy more 
than fifteen times the electric energy. 

However, results obtained from substituting 
quantum values into classical formulae do not 
mean very much. Take the example of the ratio 
of spin M to magnetic moment yz which according 
to quantum theory and experience is u/M = e/mc 
corresponding to a g factor 2. In the classical 
field theory the spin is due to the radial electric 
field of the pole ¢ and the longitudinal field of 
the dipole » which together produce angular 
momentum about the dipole axis, proportional 
to the product eu. In case of a point pole and 
point dipole with surrounding Coulomb fields, 
becomes infinite. If instead we take a “‘model”’ 
in which the field is present outside the radius 
1/k only, M becomes 2euk/3c. In our invariant 
field theory in which Maxwell and meson terms 
are subtracted the spin happens to be zero. It is 
not surprising that the classical theory cannot 
yield the correct spin, because the angular 
moment rests on the exact knowledge of per- 




































184 


pendicular components of E and H simultane- 
ously which cannot be measured without uncer- 
tainty. The actual value M=$h is just halfway 
between the field value M=0 and the ‘‘normal”’ 
value M=h. 

We learn, however, from these considerations 
that the ratio of magnetic moment to spin and 
also to magnetic energy varies widely with the 
special theory under consideration. The ‘‘model”’ 
with fields outside 1/k yields 


Wa=sek, Wireg=teh?, M=euk/c, 
1/k=23.3e/2mce. 


Quantum theory and experience, however, yield 
(see Part IV): 


Wa=4Ek, Woe=0, M=}h, 
1/k=e/2mc. (33) 


The result that the magnetic energy and mag- 
netic mass are zero in spite of the presence of a 
magnetic moment yu is suggested already by the 
simple formula n= eh/mc; if m in this formula 
should depend on yu implicitly, the ratio between 
electric and magnetic energy would have com- 
plicated values [like 16.47 (see above) ]. Also, 
the radius 1/k would be much too large as soon 
as the magnetic energy plays any considerable 
part in the self-energy [even with so small a 
factor as $ (see above) |. Already the discussion 
of various classical models and field theories 
shows that it is well possible already on a clas- 
sical basis to construct a dipole moment whose 
surrounding magnetic energy happens to vanish. 


ALFRED LANDE AND LLEWELLYN H. THOMAS 








The satisfactory result of Part 1V that the elec. 
tron does not have magnetic self-energy arising 
from its magnetic moment, and that its radiys 
is the electrostatic radius 1/k=€/2mc*, is ex. 
plained in the following way. In the first place, 
the only sources of the field according to the 
field equations (25), (25’) are the electric 
charges, and there are no retarded potentials 
from magnetic poles or dipoles. Turning to the 
Fourier representation of Part II, Section 7, the 
Hamiltonian consists of terms referring to the 
particles alone, to the field alone, and perturba. 
tion terms. In our unitary theory the terms 
referring to the particles alone vanish [first two 
terms in II, (26) ]. The part of the ‘‘mutual” 
energy which is proportional to ¢ may be inter. 
preted as electromagnetic “self-energy.” It arises 
from scalar and vector potentials, but for a 
particle at rest the contribution of the scalar 
potential [last term in first line of II, (26) ] is of 
opposite sign and half as large as the contribution 
of the vector potential [second line in I, (26); 
the resultant energy is —}ke+ke= ke. If the 
particle is in an external field one also obtains 
mutual energy terms proportional to e, namely, 
for a particle at rest e-(V’—V”’)ext and also, if 
Dirac’s quantum method is applied, the scalar 
product (u-H’ —H’) where u stands for ¢h/2me 
as though the electron had a dipole moment y, 
There are no terms, however, proportional to 
which could be interpreted as magnetic self- 
energy. In spite of displaying spin and magnetic 
moment y, the electron at rest has only electro- 
static self-energy, and the mass is m=}&k/c. 
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The accommodation coefficient of helium on platinum has been measured at 77, 193, 273, and 


373°K; and the values obtained at these temperatures for a clean platinum surface are 0.090, 
0.043, 0.071, and 0.072 (+0.004), respectively. When the platinum is allowed to remain in 
contact with the helium for twenty-four hours or more, the accommodation coefficient values 
increase to 0.43, 0.071, 0.170, and 0.170 (in the same temperature order). Little adsorption 
occurs at 193°K, but strong van der Waals adsorption appears at 77°K and activated adsorption 


at the higher temperatures. 











HE accommodation coefficient a of a gas 
striking a solid surface is defined by the 


following equation : 
a=Q'/Q, 


where Q’ is the actual average energy exchange 
and Q is the energy exchange that would occur if 
the gas molecules, on the average, attained on 
collision the temperature of the solid surface. 
Since the accommodation coefficient is not inde- 
pendent of temperature, closely approaching 
solid and gas temperatures are understood. 

The differential equation for the temperature 
distribution along a wire carrying a current and 
losing heat to its surroundings by radiation, con- 
duction through the ends of the wire, and 
conduction to gas molecules impinging on its 
surface is: 

d*t/dx* = At—B,'? (1) 


where 


A=(1/KarL)8xeorT*L 
+2(2m)!aprL(R/MT)'—L?RocL J (2) 


and 
B=[°Ro/Krr’*L (3) 


are constants. The quantities involved are defined 
as follows: 
J—current passing through wire, 
Ro—tesistance of wire at bath temperature, 
e—emissivity of Pt at bath temperature, 


_*This paper represents the publication in part of a 
dissertation submitted to the faculty of Bryn Mawr College 
in partial fulfillment of the requirements for the degree of 
Doctor of Philosophy. 

t Now with Institutum Divi Thomae, Cincinnati, Ohio. 
‘J. K. Roberts, Proc. Roy. Soc. A135, 192 (1932). 
*B. Raines, Phys. Rev. 56, 691 (1939). 
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K-—thermal conductivity of Pt at bath temper- 
ature, 
o—Stefan-Boltzmann 
watt cm~? deg.~‘, 
R—gas constant—8.314 joules °C-", 
c—temperature coefficient of resistivity, 
r—radius of wire, 
L—length of wire, 
7—temperature of surrounding bath, 
M—molecular weight of gas, 
p—pressure of gas, 
a—accommodation coefficient. 
The solution of Eq. (1) for A>0 is: 


B/A 


evAL —e-VAL 


constant—5.722 x 10-" 


XC —(1—e%44)e%47+-(1—e44)e-“47]+B/A (4) 
and the average temperature of the wire is: 
t=B/A[f(A)—1], (S) 


4—2(eV4b+¢-v4L) 
VAL(ev44 —e-¥44) 


where 





f(A) = (6) 


Upon equating the heat supplied to the wire to 
the heat losses from the wire, the following equa- 
tion results: 








f(A) 8xeorT*Lm 
i—_ 
2(2r)taprL(R/MT)'m 
+ Rec , 
where 
m= (R—Ro)/Q=Roct/Q (8) 


















is a constant for a given bath temperature 7, 
provided the power input is not too great. 

Following the method of Cox,* one observes 
that the right-hand side of Eq. (7) is independent 
of the temperature rise; and, therefore, for zero 
power input, 


f(Ao) = —14+82e0rT*Lm/ Roc 
+2(2e)laprL(R/MT)'m/Roc. (9) 


Also, for zero power input, Eq. (2) for A becomes 
8reorT*L — (2r)!2aprL(R/ MT)! 
7 Krr-L 





10) 


And it can be seen from Eq. (9) that 
— f(Ao) =1—(m/Roc)Krr-LAo, (11) 


which defines a series of straight lines depending 
on the parameter, 


N=(m/Ryc)Krr*L. (12) 


At any fixed bath temperature K and c¢ are 
known and m and R, are determined from the 
experimental data. Hence N is determined. Then, 
from the graphical simultaneous solution of (6) 
and (11), Ao can be determined. When applied to 
a run in vacuum (p<1X10~-* mm Hg) the second 
term in Eq. (10) falls out, and one can solve 
directly for the emissivity e. For a run in the 
presence of gas (1X10-*<p<1X10- mm Hg) a 
can be evaluated from Ao, once the emissivity is 
determined by a run in vacuum at the same bath 
temperature. This method of calculating at zero 
power input has increased the precision and 
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Fic. 1. Emissivity and resistance of platinum 
filament vs. temperature. 


*M. Cox, Phys. Rev. 64, 241 (1943). 


40. Reynolds, Phil. Trans. 170, 727 (1879). 
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materially reduced the work of determining the 
accommodation coefficient from the experimental] 
data. 


EXPERIMENTAL PROCEDURE 


The vacuum system was a standard gas circy.- 
lating system similar to those used by other 
workers in the field with the exception that all 
stopcocks directly in the path of the helium were 
replaced by mercury cut-offs. This was done to 
prevent contamination of the gas by stopcock 
grease. The helium was purified by high voltage 
discharge in a misch metal bottle. Pressures jn 
the experimental tube were measured by a 
McLeod gauge connecting directly to the tube 
through a necessary mercury trap. The pressures 
read were corrected for the difference in temper. 
ature between the tube and the McLeod gauge." 
Stepwise variations of pressure in the experi- 
mental tube were obtained by use of a mercury 
control. This consisted of two concentric glass 
tubes, the outer one connected to a lift containing 
mercury, the inner one containing a series of 
small holes of leak dimensions growing larger 
toward the bottom. 

The experimental tube was somewhat pot- 
bellied in shape, being 33 cm wide at its widest 
point slightly below the center. This was to allow 
sideward sag of the platinum wire when it was 
flashed. The platinum (obtained from Baker and 
Company and of 99.99 percent purity) was in the 
form of a filament (14.23 cm long at room temper- 
ature, 2.6510-* cm in radius) supported verti- 
cally by tungsten leads without the use of springs. 
Since the flashing temperature of platinum (about 
1200°C) is low, springs were unnecessary ; and 
considerably greater accuracy can be obtained by 
their elimination. 

The potential drop across the tube and across 
a ten-ohm standard resistance in series with it 
was measured with a Leeds and Northrup Type 
K potentiometer. Readings were taken with the 
current flowing first in one and then in the re- 
verse direction to eliminate errors due to thermal 
e.m.f.’s. From these readings the resistance of the 
tube and its power input were determined. 

The temperature of the medium surrounding 
the platinum was set by the use of fixed tempera- 
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TABLE I, Data under vacuum conditions. 





HELIUM 








Temperature “K 
Resistance Ro ohms 

Resistivity p ohm cm X10 : 
Thermal conductivity K watts/cm °C 


Temp. coeff. resistivity c=1/RodR/dT °CXx 10 23, 


Emissivity ¢ 


4 . 3.2 3.1 

2562 4.4172 6.4475 8.9620+0.0001 
95 6.85 10.0 13.9 +0.1 
796 0.679 0.695 0.715 +0.005 
1 5.70 3.90 2.81 +0.01 
009 0.024 0.035 0.050 +0.001 











ture thermostating baths. Ice, solid carbon di- 
oxide and acetone, liquid air, and boiling water 
were used. 


VACUUM DATA 


Initially the resistance of the wire at 0°C 
was 6.4780 ohms, and the emissivity was equal 
to 0.057 (a thermal conductivity of 0.695 
watt/cm °C is assumed).® The platinum was 
flashed first at 1000°C, then at 1200°C. At inter- 
vals the tube was thermostated and readings of 
Ry, e, and a were taken. After 100 hours of 
flashing, the ratio m/R» had become constant 
within its uncertainty, but the resistance at 0°C 
continued to decrease slowly ; and it did not reach 
stability until after a thousand hours of flashing 
had elapsed. 

Once stability had been attained at 0°C, runs 
were made at other temperatures. Reference to 
the equations will show that at zero gas pressure 
both the emissivity and the thermal conductivity 
of platinum enter. A search of the literature indi- 
cated some confusion in the values of the thermal 
conductivity of platinum. The best values ap- 
peared to be those by Holm,’ whose work was 
done on platinum of 99.95 percent purity, and 
whose results are given by: 


K=0.699[1+2.83 
X 10-4(T — 292.7) ] watt/cm°K. 


This formula was determined for temperatures 
above 0°C, so that its use at —80°C without 
substantiating data is to be questioned. But use 
of the calculated value (0.679) at this tempera- 
ture leads to a value of the emissivity which falls 
on a smooth curve between the higher tempera- 
ture values and the known value of zero at 0°K. 
The emissivity-temperature relationship is shown 
in Fig. 1. At 77.4°K the radiation loss (which is 


*R. Holm and R. Stormer, Wiss. Veroff. a. d. Siemens- 
Konzern 9, 2, 312 (1930). 





proportional to 7*) has become so low that it can 
be neglected, and the thermal conductivity is the 
only unknown. Unfortunately there is a limiting 
value of N for very small A» where the two ex- 
pressions of —f(Ao) are tangent over a con- 
siderable range, and there is quite a lack of 
precision. This value of N equals L?/2. It is 
necessary, therefore, to use slightly different 
methods of computation at this temperature. 
Two were tried—first, using a non-zero power 
input Q and the corresponding average tempera- 
ture t; second, expanding —f(Ao) in terms of 
power series in Ao and neglecting all powers of Ao 
higher than the first. The two methods gave the 
same value of K. 

Figure 1 also includes a graph of the variation 
of the resistance of the platinum wire with abso- 
lute temperature (over the range 200—400°K it is 
close to a straight line). The vacuum data are 
summarized in Table I. 


GAS DATA 


Workers'** in the field have long ago dis- 
covered that the accommodation coefficient is a 
function of the surface of the metal. In order to 
have meaning @ must be measured on the metal 
itself: That is, the surface must be free of ad- 
sorbed or combined films. A clean surface can be 
obtained by flashing the metal at high tempera- 
tures in vacuum. But, unfortunately, as soon as 
the metal is allowed to cool in the presence of a 
gas, the surface becomes contaminated. Whether 
the contamination is due to adsorption of the 
impinging molecules of the pure gas or to the 
adsorption of minute impurities from the gas, the 
effect is practically the same. The accommodation 
coefficient will increase in value as time elapses. 
It is necessary then to extend the data back to 
the time of zero contamination. Two methods 
have been used to obtain this. The first, initiated 


*W.C. Michels, Phys. Rev. 40, 472 (1932). 
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Fic. 2. Accommodation coefficient of helium on platinum 
vs. time after gas admission, 0°C, for clean and dirty 
platinum, pure and impure helium. 


by Roberts,' 7 * consists of flashing the metal in 
vacuum, letting in the gas (metal still flashing) 
until equilibrium pressure is obtained, then sud- 
denly cutting the flashing current off and taking 
readings of resistance against power input from 
that moment on. The second, used by Mann and 
Newall® and Raines,’ consists of flashing the wire 
in vacuum until a clean surface is obtained, then 
letting in the gas and taking resistance-power 
readings from the time of entrance of the gas. 
Both methods were tried in this work under com- 
parable conditions (hours aging of wire, gas 
purity, and pressure, etc.). It was noticed that at 
equal times after zero (10 min. or more) the first 
method invariably gave a higher accommodation 
coefficient than the second. On runs by the first 
method alone, a» increased whenever gas pressure 
was increased or the platinum was allowed to 
flash for longer times in the presence of the gas. 
There seemed to be only one reasonable expla- 
nation of this : namely, that residual impurities in 
the helium (despite the fact that the gas ap- 
peared pure viewed through a pocket spectro- 
scope) were adsorbed upon the platinum when it 
was flashed in their presence, and, therefore, the 
first method gave a value of a» characteristic of an 
already contaminated surface. For this reason it 
was discarded, and the second method used 
throughout the remainder of the work. 
Plotted in Fig. 2 are three curves representing 
the accommodation coefficient against time after 


7 J. K. Roberts, Proc. Roy. Soc. A129, 146 (1930). 
8 J. K. Roberts, Proc. Roy. Soc. A142, 518 (1933). 
. B. Mann and W. C. Newall, Proc. Roy. Soc. A158, 
397 (1937). 
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the entrance of the gas for three runs at 0°C 
Curve A was done with spectroscopically pure 
helium on platinum that had been flashed fo, 
many hours in vacuum before the run was takep 
and was presumably clean of adsorbed impuri. 
ties. Curve B was taken with a clean wire byt 
impure helium. Curve C was done in the following 
manner. The platinum was allowed to remain jp 
the presence of the helium until a stable accom. 
modation coefficient a, = 0.166 was reached. The 
tube was then evacuated and a resistance-power 
measurement determined for p<1X10-* mm Hg, 
Then the gas was allowed to re-enter (no flashing 
had been done) and resistance-power measure. 
ments taken against time. Curve C shows the 
resulting accommodation coefficient. It is notige. 
able that it immediately rises from a zero value to 
thé previously determined value of a,. But mere 
evacuation will not remove a film from a metal, 
Therefore the plotted values of a@ in the first few 
minutes after the gas entrance are not true 
measures of the accommodation coefficient. 
During this time the wire is cooling. from its 
equilibrium temperature in vacuum to a lower 
temperature determined by ap, and an unsteady 
state exists. The equations for the determination 
of a assume a steady state (no time term), so that 
the apparent accommodation coefficient de- 
termined by them is not true until the tempera- 
ture of the wire has ceased to lag behind the value 
determined by the energy losses. 

Reference to curve B shows, however, that the 
fastest change in the accommodation coefficient 
is occurring in the early time intervals when the 
cooling lag is present. For the above reason it was 
decided to use logarithmic extrapolation of the 
accommodation coefficient to zero time, following 
the procedure suggested by Michels.'® That this 
extrapolation is not always too good is only to be 


TABLE II. Accommodation coefficient of He on Pt at 
zero time (ao), some twenty minutes later (a,,), and after 
twenty-four hours (a..). 











Temp. 
= ao eq be 
77.4 0.090 0.42 0.43 
193.2 0.043 0.045 0.071 
273.2 0.071 0.080-0.113 0.170 
373.1 0.072 +0.004 0.110 0.170 








1° W.C. Michels, Phys. Rev. 52, 1067 (1937). 
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expected. The presence of more than one im- 
purity with different rates of film build-up is 
logical; an actual saturation effect (that is— 
entrance of the impurity within the body of the 
platinum) is possible. 

The accommodation coefficient data are sum- 
marized in Table II for the four different bath 
temperatures used. By do is meant the accommo- 
dation coefficient at zero time, and the values of 
it given in the table were obtained by logarithmic 
extrapolation except at 193.2°K where no ad- 
sorption occurred until many hours had passed. 
By @-q is meant the value of the accommodation 
coeficient obtained some ten to twenty minutes 
after the entrance of the gas. At this time the 
value of the accommodation coefficient is changing 
very slowly. By a. is meant the accommodation 
coefficient 24 hours or more after the gas has been 
let in. Figure 3 shows representative runs at each 
temperature for thirty minutes after gas admis- 
sion. Although the curves have been shown 
extrapolated to zero time the values of accommo- 
dation coefficient obtained by such extrapolation 
mean little because of the cooling lag of the 
platinum present in the first few minutes, and the 
actual values of ad» were obtained by linear 
logarithmic extrapolation. 

There are some points of interest in the ac- 
commodation coefficient data obtained. The 
values of ao at 273°K and 373°kK, i.e., 0.071 and 
0.072, respectively, are in fair agreement with the 
classical value'® of 0.078 determined from 


4mM/(m+M)?, 


where m and M are the atomic weights of plati- 
num and helium. The low value of a, (0.170) at 


TABLE III. Accommodation coefficients of 
He on Ni, W, and Pt. 








Accommodation coefficients found by other workers 

















He on Ni? He on We He on Pt 

Temp. 

°K ao Gx ao Qe ao ax 

77.4} 0.048 0.410 | 0.025 — 0.04* 
195 0.060 0.423 | 0.046 — — — 
273 | 0.071 0.360 | 0.057 0.19 | 0.03T — 

(295°K) 

373 | 0.077 0.343 — — 0.07 0.25 











* See reference 1. 

° See reference 11. 

* Bath temperature 80°K, filament temperatures 100°C and up. 

t Bath at room temperature, filament between 100°C and 1000°C. 
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Fic. 3. Accommodation coefficient of helium on platinum 
vs. time after gas admission for various bath temperatures. 


these temperatures indicates a tightly bound 
adsorbed film on the platinum and, therefore, 
activated adsorption of minute impurities present 
in the helium. At 77.4°K the high value of 
a, (0.43) combined with the fact that the best 
straight line logarithmic extrapolation was found 
at this temperature strongly indicates van der 
Waals adsorption of helium. 

In Table III are presented values of the ac- 
commodation coefficient of helium on several 
metals found by other authors. Comparison of 
the values of a determined in this work with those 
determined by Mann" (He on Pt) are difficult 
because at no time in the latter’s work did the 
filament temperatures run below 100°C. In the 
one comparable case (both bath and filament 
temperatures at 100°C) agreement is obtained for 
do (0.07 against 0.072). There seems to be little 
correlation between the present work and the 
other two sets of data (He on Ni? and W') shown 
in Table III. In the case of the nickel data which 
were determined by experimental means very 
similar to this work no attempt was made to 
allow for a possible lag in cooling of the nickel 
when the gas was admitted. None of the three 
sets of data presented in Table III shows the 
characteristics of the helium on platinum ac- 
commodation coefficient data discovered in this 
work, i.e., a region of strong van der Waals 
adsorption at very low temperatures, followed by 
a region of no adsorption and finally an extended 
region of activated adsorption as the temperature 
is increased. 


"'W. B. Mann, Proc. Roy. Soc. Al46, 776 (1934). 
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On Supersonic Waves in Cylindrical Tubes and the Theory of the 
Acoustic Interferometer 


P. E. KRASNOOSHKIN 
Laboratory of Oscillations, Physical Research Institute, Moscow State University, U.S.S.R. 


(Received May 6, 1943) 


A theory of the Pierce acoustic interferometer is given which takes into account the non- 
uniformity of the acoustic field. This leads to an interpretation of the satellite peaks observed 
in the experiments of Pielemeier and others, and to an improved formula for the determination 
of the absorption coefficient from the observational data. 





1. INTRODUCTION 


HE successful application of the Pierce 
acoustic interferometer to the determina- 
tion of supersonic sound velocities is well known. 
Many attempts have been made to apply it to the 
measurements of absorption coefficients as well, 
but the accuracy obtained has appeared ex- 
ceedingly low. L. Mandelstam has concluded that 
the most likely cause for this is the non-uni- 
formity of the sound field, for which no correction 
has been made. Ordinary methods of reduction of 
the experimental data have assumed that the 
quartz plate used as a radiator vibrates like a 
piston, giving rise to a purely plane wave. How- 
ever, real plates do not come up to this ideal, 
usually exhibiting a rather irregular form of 
vibration. 

The present paper deals with the theory of the 
interferometer, taking into account this more 
complex pattern of the acoustic field, and de- 
scribes a new method of reduction of the experi- 
mental data which is believed to be more accu- 
rate and general in application than those 
previously employed. 

Before attacking the general theory, it may be 
desirable to give a brief description of the 
principle of the interferometer and of the manner 
in which it is employed. A cylindrical metal tube 
of arbitrary sectional form is closed at one end by 
a quartz plate driven by an electrical sine wave 
generator, and at the other by a movable piston 
used as a reflector. The whole volume of gas in 
the tube is set into vibration, in virtue of which 
there is an acoustic reaction on the sound source. 
The observational result is a reaction curve in 
which some quantity measuring the acoustic 

output of the source is plotted as a function of the 


distance between the emitter and reflector plates, 
The reduction of the observations consists of the 
deduction of the ultrasonic velocity and the 
absorption coefficient from the form of this 
curve. In order to study this latter problem we 
must find an analytical expression for the acoustic 
reaction on the emitter in terms of the velocity of 
sound on the one hand, and in terms of parame. 
ters describing the form of the vessel and the 
character of the vibration of the quartz plate on 
the other. 


2. GENERAL THEORY 


A rigorous solution of this theoretical problem 
would be extremely difficult. In our analysis we 
shall consider the various types of waves which 
may be generated in such a system, and their 
effect on the output of the source. This mathe- 
matical theory is similar in principle to the 
theory of electrical wave guides. 

If, ata given frequency, the absorption may be 
considered as small and the transverse dimen- 
sions of the cylinder as large, we may follow 
Kirchhoff! in neglecting viscosity and heat losses 
on the surface of the cylinder, and consider it as 
infinitely rigid, perfectly smooth, and absolutely 
thermonegative. By virtue of these conditions, 
we assume the existence of a velocity potential 
satisfying the equation 


Vee tkio=0 (1) 


throughout the interior of the cylinder. Here 
ky=Bo—tao is a complex wave factor for the 
unconfined medium (89=w/vo), and the time 
factor e** is understood. 

With the z axis parallel to the generators of the 


1G, Kirchhoff, Pogg. Ann. 132, 177 (1868). 
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cylinder, the reflecting surface in the plane =0, 
and the radiating surface in the plane z=L, the 
boundary conditions take the form 


dy/an=0, on the side walls of thecylinder, (2) 
ay/dz=0, in the plane z=0, (3) 
ag/ds= F(x, y), in the plane z=L, (4) 


where F(x, y) gives the normal component of the 
velocity on the surface of the radiator. 

Following Rayleigh’s method? we look for 
particular solutions of Eq. (1) in the form of 
progressive plane waves propagated in the z 
direction, by setting 


g=V(x, ye=m™, (5) 
where y satisfies the condition 
ay day 
—+—-+ (ko? —m?)y=0 (6) 
dx? dy? 


with f=ko?—m?*, and where it will be assumed 
that f is real. 

Supposing this equation solved, with boundary 
condition (2), we find the types of waves which 
may be propagated along the cylinder without 
change of type, corresponding to a discrete set of 
values f;, 7=0, 1, 2, --- for f, with associated 
values for the wave factor m;=8;—ia;. Under the 
conditions assumed, the different solutions of 
Eq. (6) will form an orthogonal set of functions 
when integrated over the cross section of the 
cylinder, and we shall suppose further that they 
have been normalized. 

The types of waves permitted vary with the 
form of the bounding cylinder, as do their phase 
velocities and attenuation coefficients. The wave 
fo=0 corresponds to an ordinary uniform plane 
wave ; the greater being the value of f;, the more 
complicated the wave pattern becomes. 

Separating the real and imaginary parts of 
m?=k,?—f;, we obtain 


28? = (Bo? — ao* — fj) 

+[ (Bo? — a0? — f)?+401078o"]#, (7) 
2a;?= — (Bo? — ao? — fj) 

+[ (80? — a0? — fj)? +-4c07Bo?}#. (8) 


*Lord Rayleigh, Theory of Sound, Vol. 2, Section 301. 
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The typical form of dependence of 8; and a; on 
the value of f; is shown in Fig. 1. The waves for 
which f; is of the order of 80? have such large 
attenuation that they practically disappear 
within a few wave-lengths, while if a»=0, waves 
for which f;2 Bo? are not propagated at all. This 
type of attenuation is not due to adsorption in 
the medium, but to interference effects of pri- 
mary waves from the source and their reflected 
waves from the sides of the cylinder. It is this 
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Fic. 1. Attenuation and phase velocity of corrugated waves 
in tubes for ao/Bo=6X 10~. 


effect which constitutes the basis for the present 
theory. 

The general solution of Eq. (1) is obtained as a 
linear superposition of all of these special waves 


e=L vi(x, »)[c;’ exp (tmz) 
+c; exp (—imgz)]. (9) 
From the form of this solution, boundary condi- 
tion (2) is satisfied automatically. Condition (3) 
will be satisfied by taking cj =c;'’ =c;. 
We assume that, whatever the form of vibra- 


tion of the source, F(x, y) of condition (4) can be 
expanded in a series of the form 


F(x, 9) = Awix, y) (10) 


from which, restoring the time factor, our general 
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Fic. 2. The elementary reaction curve for air at a frequency of 100,000 cycles. 


solution satisfying all conditions turns out to be 


A Vi 


1m; 


exp (imjz)-+exp (—im;z) 





y= Rj e*! > 





» (11) 

exp (im;L) —exp ee 

where the symbols ® and g will be used to in- 

dicate the real and imaginary parts, respectively. 
The acoustic output of the radiator is 


dw d¢g d¢g 
f°) (2) os 
in which p is the density of the medium. The 
integral is to be extended over the surface of the 
radiator. 
By substitution of Eq. (11) in (12), taking 
advantage of the orthonormal character of the 


y’s, and performing a time average in the usual 
fashion, we get 


dw 1 A; 
ae 
dt Av 2 B?+a; 
X[B;J(cot mL) +a;R(cot m;L)]. (13) 
To a sufficient approximation, the terms con- 
taining the factors a;R(cot m;L) can be neg- 


lected, as they are small for both small and large 
values of a;, since for large values of a; 


R(cot m;L) ~exp (—2a;L). 


In the remaining expression we need consider 
only the terms for which |8;|>|a;| in which we 
can put 


wB;/ (BP? +a;*?) ~w/B,;=0;. 


If the variation in the output of the sound 
generator over a resonance curve is small, we may 
assume that the experimentally measured plate 
current of the power tube which drives the 
oscillator is a linear function of the acoustic 
power output, and we can take as a measure of 
the observed curve the remaining quantity 


dw 
(—) 21(L)=¥ pv,Afg(cot mL). (14) 
dtl w 


3. DISCUSSION 


The simplest case is that in which the radiator 
vibrates like a piston, as is usually assumed, so 
that only Ao differs from zero. This reduces I(L) 
to the form 


I(L) = pvpAo*[1+exp (—2acL)-cos 2BoL] (15) 


for large values of L. The reaction curve for the 
ordinary theory is shown in Fig. 2. According to 
Eq. (14) the complete reaction curve is made up 
of a superposition of curves of this type.*~* 

In order to explain the usual form of observed 


3 W. H. Pielemeier, Phys. Rev. 34, 1187 (1929). 
‘J. S. Hubbard, Phys. Rev. 38, 1011 (1931). 
5 L. Beliavskaya, Bull. Acad. Sci. U.S.S.R., No. 7 (1935). 
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reaction curve, it is necessary to assume that in 
most cases the acoustic field set up consists of a 
group of waves forming a narrow wave packet. 
The presence of absorption diminishes the re- 
solving power of the interferometer, so that the 
fine structure remains unresolved. 

In particular cases, in which there are two 
dominant groups of waves, satellites may be ob- 
served, which gives us an explanation of the 
subsidiary peaks which have been observed by 
many authors.* *-* Figure 3 represents such a 
curve taken from the work of Pielemeier.* The 
measurement of the intervals indicates that the 
major peaks correspond to waves of high index 
f;, while the minor peaks correspond to the nearly 
plane wave. This is peculiar to Y-cut crystals ; for 
X-cut crystals the approximately uniform plane 
wave predominates. 

Since only the peaks without satellites are used 
in the experimental work on absorption coeff- 
cients, we confine our attention to them. Ac- 
cording to our theory, we consider them to be 
made up of a group of waves forming a single 
wave packet. We shall now show that the ap- 
parent absorption coefficient of such a wave 


(1931) Klein and W. D. Hershberger, Phys. Rev. 37, 760 
7R. Alleman, Phys. Rev. 55, 87 (1939). 
*W. H. Pielemeier, J. Acous. Soc. Am. 9, 600 (1938). 


CORRUGATED WAVES 


Fic. 3. The satellites after Pielemeier. 












IN TUBES 


3gth 
Peak 

















packet arises from two independent causes: 
(a) the true absorption of the medium, and (b) an 
apparent absorption arising from the mutual 
interference effects of the waves comprising the 
packet. The separation of the two effects will lead 
us to our improved formula for the reduction of 
the experimental data. 

For large values of aol, when the dominating 
part of the acoustic reaction on the source may 
be attributed to the first reflected wave, ex- 
pression (14) takes the form 


I(L)= 6{1+exp (—2a;L) cos (26;L)] (16) 
with 
0;= pv;A ;. 


Replacing the line spectrum by a continuous 
wave group, we get 


0 


I1(L)= | 0(8)[1+e-2" cos (28L) ds. 


Bo 


If we transfer the origin to the center of the 
packet 8., and assume the whole spectrum to be 
confined within the interval 8.+), we get 


I(L)2=C'(L) 
Xexp (—2a-L) cos (28.L+¢)+const., (17) 






















where m 
C’(L) cos o-f 6,(é) -cos 2Lé-dé, 
a 
(18) 
+b 
C’(L) sin o-f 6.(€)-sin 2LE-dé, 
— 
6.(E) = $1 0(€) +0(—£) J, (19) 


Oa(&) = 3L0(E) —0( —é) ]. 


These equations may be inverted in the form 


6,(€) =cos ef cw -cos 2§L-dL, 

“See 

3 (20) 
@.(¢)=sin ¢- f C'(L)-sin 2¢L-dL. 

0 


It is evident that the relation between the form 
of the spectrum 6(€) and that of the reaction 
curve envelope C’(L) is similar to that con- 
necting the optical spectrum and the visibility 
curve of the Michelson interferometer. 

Since a large variety of forms of the function 
6(€) may produce effectively the same enveloping 
function C’(L), it is not surprising that observed 
reaction curves may be approximated rather well 
by an over-all exponential decay curve of the 
form e~*24, 


4. REDUCTION OF THE EXPERIMENTAL DATA 


It is customary, in the reduction of experi- 
mental observations, to choose those reaction 
curves which may be represented with sufficient 
accuracy for large L by an expression of the form 


exp (—2a.L)-cos (28.L+¢). (21) 


By comparison with the theoretical curve (15) 
for a uniform plane wave, it is normally assumed 
that (a) the experimental decay coefficient a, is 
equal to ap for the uniform wave, and (b) 8,=8o, 
so that the wave-length of the uniform wave is 
equal to twice the distance between adjacent 
maxima. 

By comparison with our more complete formula 
(17), we see that the dependence on L is de- 
termined in part also by the variation of the 
envelope function C’(L). If we suppose that the 
over-all variation of the reaction curve is expo- 
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nential in character, we must put 


C'(L) = Co-e-*% (22) 


so that the apparent decay constant becomes 
ag ao+a. (23) 


The absorption coefficient a, arises from direct 
absorption in the medium, while the absorptidp 
represented by the coefficient a is due to inter. 
ference effects within the wave packet, the 
resultant absorption being due to the combina. 
tion of the two effects. 

Substituting (22) in (20) we find 





6,(é) =Ccos ¢° 
a?+é 
so that 
a=Aé 


is the half-width of the symmetrical part of the 
packet. 

Now we can express Aé in terms of the width 
Af, for according to (7) and (8) 


2B.AE=Af 
so that to a first approximation 
a=Af/8B.=Af/8Bo. (24) 


The value of Af is characteristic of the form of the 
vibrations of the source and the geometry of the 
interferometer, but does not depend on the 
properties of the medium in it. 

According to Eq. (8), if f./82<4, we have 


aa 1 + (f-/2Bo?) 1, (25) 


so that, if the width of the symmetrical part of 
the packet Af and the position of its center f, are 
known, the coefficient of attenuation ap referred 
to the supposed homogeneous plane wave may be 
deduced from the apparent absorption coefficient 
by means of the formula 


ae— (Af/8Bo) 
ao= . 
"1+ (fe/2B0?) 
while the corresponding wave-length is given by 
the formula 
Ao=A-/L1+(f-/2B0*) ] (27) 


in terms of the apparent wave-length X,. 


(26) 
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Since many authors reduce their results to a 
calculation of aoc, we have the corresponding 


formula , 
ii iin Lae — (Af/8Bo) ]-d. (28) 
(1+ (fe/2B0") J 


As far as can be judged, the corrections in- 
volving Af and f. are appreciable, and with an 
ideal uniform plane wave the ultrasonic inter- 
ferometer would give a much less rapid attenua- 
tion than that usually measured. 

The determination of Af and f, for a given 
apparatus may be carried out by observation of 





the reaction curve of a gas with known values of 
vo and ap. The coefficient a obtained by the 
method of Lubny-Herzyk® may be used in (24) to 
determine Af. With Af known, f. may be de- 
termined by measuring a, and then using (25) 
with Ape =ap—a. 

In conclusion I should like to express my deep 
gratitude to L. Mandelstam, under whose guid- 
ance this work has been done, and to K. 
Theordortchik for his exceptional attention to 
the present investigation. 


*E,. Pumper, J. Phys. U.S.S.R. 1, No. 5-6 (1939). 
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The intensities of the x-ray reflections from a given crystal of a known substance depend on 


the vector distances between the atoms and not directly on the coordinates of the atoms. The 
problem of uniqueness in the x-ray analysis of a crystal structure thus depends on the uniqueness 
of the determination of the arrangement of a periodic set of points in space by its vector distance 
set. In this paper a large number of cases is presented in which two, and sometimes three or four, 
non-congruent sets of points are homometric, i.e., have the same vector distance set. Although 
the investigation is largely based on a discussion of one-dimensional cyclotomic sets, i.e., those in 
which the points divide the period on the line into rational fractions, it is shown that there are 
many families of non-cyclotomic pairs and multiplets and that each of these families has its 
counterpart in two and three dimensions. The significance of these results for practical crystal 


analysis is discussed. 


INTRODUCTION 


T is well known that the F factors of x-ray 

analysis are completely determined in magni- 
tude and in phase by the coordinates of the atoms 
in the crystal and by the f factors of the indi- 
vidual atoms. The absolute values of the F 
factors are, however, given by 


|F(h)|*=EE f.feexp 2wih-(ee—n,), (1) 


in which f, and r,, respectively, are the f factor 
and the coordinate vector for the sth atom and 


h= hb. (1a) 
* A preliminary report of this work was presented at the 


Ann Arbor Meeting of the American Society for X-Ray 
and Electron Diffraction in June, 1943. 









In this latter equation the quantities h; are the 
Miller indices for the plane in question and the 
vectors b; are the reciprocal lattice vectors. The 
intensities in the x-ray diffraction pattern thus 
depend only on the vector distances between 
atoms and on the f factors, since in Eq. (1) the 
coordinate vectors r, do not enter alone, but only 
in the form of the differences r;—r,. 

We shall not concern ourselves here with the 
problem of determining the distances from the 
values of | F(h,;)|* since this has been discussed 
in detail elsewhere.! We shall investigate the 
problem of the determination of the coordinate 

1A. L. Patterson, Zeits. f. Krist. 90, 517 (1935); D. 
Harker, J. Chem. Phys. 4, 381 (1936); M. Avrami, Zeits. 


f. Krist. 100, 381 (1939); S. H. Yu, Nature 149, 638 (1942) 
and 150, 151 (1942). 
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Fic. 1, Sample periodic set on a line, and the representation 
of its distance set on a circle. 


vectors r, from the values of the vector distances 
f.—f,, since it is apparent that a knowledge of 
these distarices is the most information which we 
can possibly obtain from the x-ray diffraction 
data alone. 

It is well known that the methods of crystal 
analysis introduce a center of symmetry and 
that two structures which can be related to one 
another by a center of symmetry give the same 
diffraction effects. Any two structures which 
are congruent in the sense that they can be 
brought into coincidence with one another by a 
translation, a rotation, a reflection, or some com- 
bination of these operations give the same x-ray 
pattern. We shall then say that a structure analy- 
sis is unique if all the structures which can be 
derived from it are congruent to one another. 

That the distance vectors do not determine all 
structures uniquely is known from the work of 
Pauling and Shappell? who found that the special 
point 24(d) of the space group 7),’—Ja3 (for the 
values +x and —<x of the parameter x) presents 
two non-congruent structures which have the 
same x-ray intensities and consequently the same 
distance vectors. It has been shown? that this 
property is also exhibited by a corresponding 
special point for the space group O,!'°—Ja3d. 
Most of the further preliminary discussion of the 
problem*‘* must now be withdrawn, however, in 
the light of the results presented here. 

We are thus concerned with the discussion of 


( " 30) Pauling and M. D. Shappell, Zeits. f. Krist. 75, 128 
1 , 

3A. L. Patterson, Nature 143, 939 (1939). 

4A. L. Patterson, Phys. Rev. 55, 682 (1939). 
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the cases of ambiguity which may arise in the 
determination of the arrangement of a set of 
points from a knowledge of the vector distances 
for these points. We shall omit discussion of the 
finite problem, and confine our attention to the 
case of an infinity of points arranged periodically, 
We first consider periodic distributions of points 
on a line; later we shall refer to problems in two 
and three dimensions. 


DISTANCE ARRAYS 


Consider a set of N line lattices of points repre. 
sented by the coordinates x;+kM (¢=1, 2, ---N) 
with & any integer. These N sets of points can be 
equally well represented by the symbols x, 
(mod M), and together form a periodic set of 
period M. The discussion of periodic sets of this 
type is facilitated by making the transformation 


0,;=22rx;/M. 


We can then represent the set on the line by 
points on the circumference of a circle, the 
angular coordinate of the ith point being 4,, 
A sample set on the line and its representation on 
a circle are shown in Fig. 1. 

It is easy to see that the distances A ;; (mod M) 
for a periodic set on the line are given by the 


TABLE I. Survey of cyclotomic sets (r points on a period 
divided into m equal parts). The number of possible 
arrangements was obtained by first considering the number 
of r partitions of m. The number of cyclic arrangements 
(clock order neglected) of each r partition was then con- 
sidered. In all, 2664 cases (4<r<mn/2, n=16) were in- 
we. 390 pairs, 7 triplets, and 3 quadruplets were 
ound. 














Possible arrangements Number of pairs 

mir=12 3 4 5 6 7 8 4 5 6 7 s 
8 a4 § & 5 4 1 iji -—- —- — = 
9 14 #7 10 10 7 4 1jo —- — — — 
10 1 5 8 16 16 16 8 sijio0 3 - —- = 
11 15 10 20 246 26 20 10j;0 0— — =— 
12 16 12 299 38 SO 38 29)1 3 15 _ =— 
13 1 6 14 35 57 76 76 £57)}1 0 1 _ — 
14 1 7 16 47 79 #126 133 126/00 6 6 48 - 
15 1 7 19 S56 111 185 232 232};0 §5§ 23 9° — 
16 1 8 21 72 147 280 375 440|2 10 316 42% 18% 











« 3 triplets are here counted as pairs. If all pairs of the triplets are 
counted, 6 pairs should be added. 

> 4 triplets are here counted as pairs. If all pairs of the triplets are 
counted, 8 pairs should be added. 

¢ 3 quadruplets are here included as far as their complementary pairs 
are concerned. If the four non-complementary pairs of each quadruplet 
are included 12 pairs should be added. 


5 Many results concerning this and related problems are 
discussed by S. Piccard: Sur les ensembles des distances des 
ensembles des points d'un espace euclidien, Mémoires de 
l'Université de Neuchatel, 13 (1939). 
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expressions 
Aijg=Xj— Xs (i,j=1, 2, °°", N). (2) 
We shall define as a simple translation a; (mod M) 


the translation which joins any point x; with the 
point next in order of numbering. It then follows 


that 
Q;=Xin1— X= Ag 41 (3) 
and hence that 
Aijg=QitGinit:+++aj1 (t<j). (4) 
All the distances A;; for a four-point set such 
as that illustrated in Fig. 1 can then be exhibited 
in the array: 


Qi:+a2 a1+42+43 


0 a, 
—a; 0 ae A2+a3 5 
—d,—d2 —d2 0 a3 ( ) 
—Q;—G2—-a3 —d2—a3 —d3 0 


in which all operations are to be carried out 
(mod M). 

We call an array of the type of (5) a distance 
array, and now list a few of its properties: 


(A) The ith row represents the coordinates of 
the points of the set referred to the ith atom as 
origin. All rows thus represent the same set of 
points. 

(B) The ith column represents the coordinates of 
the set obtained from that of the ith row by an 
inversion in the ith point. All columns represent 
the same set of points and one which is the in- 
verse of that represented by the rows. 

(C) Every entry in the (¢+1)th row is obtained 
from the corresponding entry in the ith row by 
subtracting aj. 

(D) Every entry in the (¢—1)th column is ob- 
tained from the corresponding entry in the ith 
column by subtracting a,. 

(E) The distance array obtained by interchang- 
ing the ith and jth row and then interchanging 
the ith and jth column will represent the original 
set with a different numbering of the points. 

(F) Any distance array is fully determined by 
any row or any column combined with a knowl- 
edge of the modulus M. 
(G) Any distance array is determined by its 
simple translations a; provided their order is 
known. 

(H) Ifa set A is represented by a given distance 
array, then any set B congruent to A can be 











dy « 


Fic. 2. A simple homometric pair. The distances are the 
same but the points are in different arrangements. 


represented either by the same array or by the 
array obtained from it by interchanging rows 
and columns. 


We shall say that two arrays are equivalent if 
they represent congruent sets. It then follows 
that any array can be transformed into any 
equivalent array by some succession of the 
operations referred to under (E) and (H). 

If we are given the distances of an array such 
as (5), but have no knowledge of their order 
within the array, it is always possible to find at 
least one arrangement of these distances in an 
array of the type (5). Routines can be established 
which will do considerable preliminary sorting. 
This sorting can be completed by trial and error 
methods which must in the end be successful 
since the numbers to be arranged came originally 
from a distance array. In general this sorting 
process will have several possible solutions, corre- 
sponding to different numberings of the set. If all 
the arrays which can be obtained from a given 
set of distances are equivalent, we can then say 
that the structure defined by these arrays is 
uniquely determined by its distances in the sense 
that we have defined uniqueness for an x-ray 
crystal analysis. That this uniqueness does not 
hold for all sets on the line is demonstrated by 
the counter-example exhibited in Fig. 2. The two 
sets shown there have the same set of distances 
but are not congruent, hence the distance array 
for this set is not unique. We say that two 
arrangements which are non-congruent but pos- 
sess the same set of distances form a homometric 
pair® of arrangements. 


CYCLOTOMIC SETS 


The discovery of the homometric pair illus- 
trated in Fig. 2 suggested that further pairs of 
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Fic. 3. The cyclotomic sets n=8, r=4, It will be noted that the third and 
fourth sets of the top line are complementary and therefore form a pair (that of 
Fig. 2). All other sets are congruent to their complements and no two have the 


same distance set. 


the same type might be found by arranging 
points on the circle at r of the m vertices of an 
inscribed regular polygon. Sets obtained in this 
manner we call cyclotomic sets. 

The results obtained in the examination of 
2664 cyclotomic sets are exhibited in Table I. 
We have found 390 homometric pairs, 7 sets of 
homometric triplets, and 3 sets of quadruplets. 
These figures are by no means final although the 
examination was quite carefully carried out. No 
theory seems to exist for the number of non- 
congruent cyclotomic sets, and it is very easy to 
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Fic. 4. Three homometric triplets from the cyclotomic 
sets: =16, r=6. The three sets in each horizontal row 
form a triplet. Note the invariant subsets of each triplet, 
and of each pair in each triplet. 











omit one or more cases despite a careful systema- 
tization of the search. 

In the course of the examination of the cyclo- 
tomic sets a number of interesting results have 
been obtained. Two of these results must be 





Fic. 5. Four homometric triplets from the cyclotomic 
sets: n=16, r=7 (full circles). Most of the distance lines 
are omitted for clarity. An invariant subset for each tri 
is outlined in full lines. Variant points are connected by 
dotted lines. The open circles will form four homometric 
triplets from the cyclotomic sets: n = 16, r=9. 
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Fic. 6. Three homometric quadruplets from the cyclotomic sets: n= 16, r=8 
(full circles). Most of the distance lines are omitted for clarity. An invariant subset 
for each quadruplet is outlined in full lines. Variant points are connected by dotted 
lines. Examination of the open circles will show that the first two and the last two 
of each quadruplet form complementary pairs. 


mentioned here, but the proofs will be left for a 


more detailed discussion elsewhere. 

If we consider a set A consisting of r points 
located on an m cyclotomy, then the set A’ con- 
sisting of the m-r points of the cyclotomy which 
are not occupied by A is called the comple- 
mentary set to A. It can be shown (i) that if 
two sets A and B form a homometric pair then 
A’ and B’ also form a homometric pair, and (ii) 
that every cyclotomic set with m even and 
r=n/2 is either congruent to its complementary 
set or it forms a homometric pair with it. 

The first of these results cuts down consider- 
ably the number of cyclotomic sets which need 


’ Fic. 7. Homometric pair (4 points) with one variable 
Tepes _ by generalization of the set of Fig. 2. 
=la< 








investigating since the homometric pairs for 
(n—r) points on an m cyclotomy can be deduced 
from those for r points. The second of these 
results enables us to determine by inspection 
many of the homometric pairs for r=n/2. The 
case of r=4 and m=8 is illustrated in Fig. 3. 
There is one homometric pair (that of Fig. 2) 
and the two members are complementary to one 
another. All other sets are congruent to their 
complements and no two have the same dis- 
tance set. 

Figures 4 and 5 exhibit the triplets found on 

= 16 for r=6 and r=7, respectively, and poe 6 
ene the quadruplets found on n=16 for r=8 
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Fic. 8. Homometric pair (m+-2 points) with one variable 
Firat a: obtained by generalization of the set of Fig. 2. 








=1, a+b= (2m); a<(2m)~.] The set is illustrated 
r m= 3 but it may be constructed on any regular polygon. 
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TABLE II. Generalized five-point pairs. Each horizontal line gives the simple translations for one pair. 
a and 6b are variable parameters. 














a1 a: a3 a as ay’ a?’ a;’ a,’ a;’ 

a b 1/2—2a—2b 2a+b 1/2—a a 1/2 b 1/2—2a—2b a+b 

a b 1/2—2a—b 2a+b 1/2—a—b a a+b 1/2—b b 1/2—2a-) 
a 1/5—a a 2/S—a 2/5 a 1/5—a 1/5 a 3/5~¢ 
a 1/5—a 1/5 1/5+a 2/5—a a 1/5 1/5—a 1/5 2/5 

a 1/5—a 1/5+a 1/5 2/5—a a 1/5 1/5 1/5—a 2/5 

a 1/6 1/6—a 1/3 1/3 a 1/6—a 1/6+a 1/3 1/3~¢ 
a 1/4 1/4-—a 3/8 1/8 a 1/4—a 1/8+a 1/8 1/2~¢ 
a 1/8 1/8—a 1/4+a4 1/2—a a 1/8—a 1/8+a 1/4—a 1/2 

a 1/4—a 1/4+a 3/8 1/8—a a 1/4 1/8—a 1/8 1/2 








—== 





The appearance of these quadruplets indicates 
‘that homometric pairs may be formed with sets 
other than the non-congruent complementary 
sets. 

In many of the multiplets so far studied there 
are subsets which remain unchanged in going 
from one member of the multiplet to another. 

‘ These invariant subsets can easily be recognized 
by the reader in the sets of Figs. 2 and 4. In 
Figs. 5 and 6, in which most of the distance 


lines have been omitted for clarity, invariant. 


subsets are outlined in full lines, while the 
variant points are connected by broken lines. 
The choice of an invariant subset in a multiplet 
is not necessarily unique. 


GENERALIZATION OF THE CYCLOTOMIC SETS 


The study of the cyclotomic sets soon made it 
apparent that these were special cases of much 
more general types of pairs. The four-point pair 
of Fig. 2 immediately permitted generalization® 
in the manner shown in Fig. 7. We were thus led 
to an infinity of cyclotomic pairs connected by a 
continuum of pairs which do not fall in the 
cyclotomic class. Another type of generalization 
based on the pair of Fig. 2 is exhibited in Fig. 8. 
It can be developed from any regular polygon 
and can thus occur for (m+2) points. Many 
other general types of pairs have been discovered 
and only a partial report can be given at the 
present time. Table II lists the five-point families 
found by generalization of the cyclotomic pairs 
of Table I. It may well be incomplete. Over 
thirty six-point families have been discovered by 
an incomplete investigation of the pairs of 
Table I. Their listing must be postponed, how- 


* This generalization was suggested by Dr. P. Erdés who 
also provided helpful discussion of several other points. 


ever, until an adequate classification system has 
been developed. The triplets and quadruplets of 
Figs. 4—6 also produce families of multiplets, but 
the detailed discussion is incomplete and must 
be withheld at the present time. 

An exhaustive study of all possible pairs and 
other multiplets seems to be of extreme difficulty, 
There is obviously a limit to the number of 
families of multiplets which can occur for a given 
number of points, but so far it has not been 
possible to set up a theory which will establish 
this limit. 


EXTENSION TO TWO AND THREE DIMENSIONS 


Little has been done so far on the extension of 
the study to two and three dimensions. It is 
easy to show that every one-dimensional family 
of pairs has its counterpart in two- and three- 
dimensional periodic sets. As a modulus M we 
can use any translational group or translational 
subgroup of any lattice in the plane or in space, 
In the construction of the distance array we use 
the vectors a, @2, &3, etc., instead of scalars and 
carry out all our calculations vectorially using 
M as modulus. Corresponding to every homo- 
metric pair which we can obtain on the line there 
will be at least one homometric pair in space. 
We are assured of this because every translation 
group in space is either a line lattice, or it con- 
tains infinitely many line lattices as subgroups. 
We can use any such line lattice as a modulus 
M’ and calculations (mod M’) will also be valid 
(mod M). The vectorial calculation (mod M’) is 
then exactly parallel to the scalar calculation 
(mod M). In the pairs (mod M) we simply re- 
place fractions of M by fractions of M’, and 
replace disposable scalar parameters by dis 
posable vector parameters. These new solutions 












































7 @® se 


Qa.ec Qwnm 





SIONS 


ion of 
It is 

















are by no means trivial since the vector param- 
eters are no longer confined to the direction of 
the line lattice M’. An example of a two-dimen- 
sional pair in which the modulus M’ is the 
diagonal line lattice of the primitive cell is 
exhibited in Fig. 9. It is based on the one- 
dimensional family of Fig. 7. 

It is obviously a restriction to use as modulus 
a translation subgroup M’ rather than the full 
translation group M. If this restriction is re- 
moved homometric pairs may appear in the plane 
or in space which have no counterpart on the 
line. We have been able to show that the two 
three-dimensional pairs already discussed** can 
be based on a two-dimensional modulus M’ which 
is a subgroup of a three-dimensional lattice M. 





CONCLUSION 


The picture presented in this paper is some- 
what pessimistic as far as x-ray analysis is con- 
cerned. It must be remembered, however, that 
in very few cases are the atoms of a crystal all 
of one kind, and it seems very probable that the 
presence of a second kind of atom will often’ 
resolve ambiguities which might occur in the 
location of the first if taken alone. A detailed 
examination must ultimately be made of all the 
cases in three dimensions in which ambiguities 
can occur for arrangements of one and of several 
kinds of atoms. 

Before a particular structure can be deter- 
mined by x-ray data alone, a complete examina- 
tion of all ambiguities which can arise in that 
structure must be available. In addition, it must 
be remembered that, owing to experimental in- 


7 Pauling and Shappell, reference 2, have shown that a 
second kind of atom in certain special positions of 7,’—Ja3 
does not resolve the ambiguity exhibited by that group. 
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accuracy, x-ray analysis does not usually give all 
the distances or information equivalent to a 
knowledge of all the distances. Such a lack of 
information can lead only to further ambiguity. 
It seems therefore that x-ray analysis must be 
used in combination with all other available 
physical, chemical, and morphological data in the 
determination of the structure of the solid state. 
That x-ray analysis can in special cases stand 
alone is undoubted,® but each of these special 
cases requires its own investigation, not only of 
the geometry of the situation in question, but of 
the accuracy of the methods in use, to avoid the 
possibility of ambiguity. 

The writer wishes to express his indebtedness 
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Fic. 9. (a) and (b). Homometric pair in two dimensions 
derived from the pair of Fig. 7. The modulus M’ is the 
diagonal of the unit cell. The lattice M is not necessarily 
orth - (c) Distance vector map for the pair of (a) 
an , 


to Mr. John C. Oxtoby of the Department of 
Mathematics of Bryn Mawr College for his 
continued interest and his invaluable aid in the 
development of the mathematical background 
for this paper. 


® Reference is made of course to the structures of many 
elements and simple compounds and to such investigations 
of more complicated compounds as those of Dickinson and 
Raymond on hexamethylenetetramine, West on potassium 
dihydrogen phosphate, and Robertson on the phthalo- 
cyanines. 
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Anomalous Rate of Nuclear Disintegration 
Effected by Cosmic Rays 


A. P. Zupanov, N. A. PERFILOV, AND M. Y. DEISENROD 
Radium Institute of the Academy of Sciences of U.S.S.R., Kazan, U.S.S.R. 
January 26, 1944 


HE disintegration of atom nuclei by cosmic rays may 
be only rarely observed in the Wilson chamber, and 
the method of thickly coated photographic plates first 
proposed and worked out! at the Radium Institute of the 
Academy of Sciences of U.S.S.R. has proved to be more 
effective. Applying this method to a series of experiments,? 
we succeeded in observing many individual cases of nuclear 
disintegration. The rate of nuclear disintegration at sea 
level is about 8X 10-* an hour per 1 cm? of the photoplate 
area and shows a rapid increase depending on altitude and 
reaching, for example, at 7000 m, an amount 50 times as 
large, as that at sea level.? This paper records some pre- 
liminary data obtained in the case of an exceedingly high 
rate of nuclear disintegration we observed late in 1942. 
The investigation was made in the following way: 
Thickly coated (70u) plates with the emulsion sensitive to 
the charged particles analogous to the emulsion E; applied 
so far' were prepared for a series of tests to be carried out 
on November 19, 1942. The plates developed three days 
after their preparation (and not radiated at all) upon being 
examined under the microscope were found to contain an 
enormous quantity of forks, which because of a great range 
of particles had necessarily to be attributed to the nuclear 
disintegration from cosmic rays. This conclusion is also 
supported by the outward appearance of tracks, for the 
emulsion E; as it was formerly reported’ possesses a 
property contributing by their appearance to the dis- 
crimination of tracks of mesotrons, protons, a-particles, 
and other particles capable of higher ionization. It was 
noted earlier‘ that the property of this emulsion allowed us 
to point out the possibility of registering mesotrons and of 
evaluating the mass of these particles. Yet in the plates 
developed on the day following their preparation, the 
number of similar disintegrations fully coincided with the 
above-mentioned “‘normal’’ rate for sea level. Further 
plate development (within different intervals up to 1} 
months) led to the observation of a greatly increased 
number of nuclear disintegrations. The examination of 
plates also showed the presence of a whole series of charac- 
teristic anomalies. Thus, for example, the former studies 
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(carried out both at our Institute and in other laboratories) 
have shown that at sea level there prevail disintegrations 
involving three heavy charged particles, the number of 
ejected particles slowly increasing with the altitude; the 
quantity of single heavy tracks is much higher than that 
of forks. In our case the greatest percent of disintegrations 
is obtained from fourfold forks, and the quantity of single 
tracks is obviously smaller than that of disintegrationss 
Side by side with single tracks and forks of a small number 
(2-5) of tracks one may observe (mostly on plates devel. 
oped in December) an increased number of disintegrations 
with the ejection of many particles and characteristic 
showers consisting of a great quantity (up to some score) 
of heavy particles that are analogous to the shower of 
about 100 particles which was registered in 1938 on the 
plates sent up in an airplane to the altitude of 9000 m.4# Ay 
illustration thereof is given in a stereomicrophotograph of 
disintegration with the ejection of 20 heavy charged par. 
ticles in Fig. 1 and an example of a shower of 50 particles 
(protons) in Fig. 2. 

It is necessary, however, to point out the great flucty. 
ations in the distribution of nuclear disintegrations on the 
plate area. For example, on some plates with the area of 
10 mm? the magnitude of the rate amounts to the value 1@ 
as high as that of the aforesaid ‘‘normal’’ value. Neverthe. 
less, the greatest rate seems to have taken place during the 
period from November 23 till December 1. But it is only 
after a detailed and painstaking investigation of the plates 
under a microscope that we may come to a final conclusion, 


20 CM AIR EQUIV. 








Fic. 1. Example of disintegration involving 20 heavy particles, 
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Fic, 2. A shower of 50 particles (protons). 
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The question becomes complicated due to the presence 
of an abnormal concentration of forks in the proximity of 
big showers—the fact noted at an earlier* period and 
graphically seen when drawing a topographic map of the 
distribution of nuclear disintegrations.* Doubtlessly the 
reason that called forth this anomalous phenomenon is a 
rough increase of the intensiveness of those components of 
cosmic radiation that are responsible for nuclear disin- 
tegrations. Considering this, it is of great interest to find 
out whether the increase is a general one and how long it 
will last (it may be the result of a cosmic catastrophe and 
thus a similar anomaly could be observed in other places 
of the globe), or whether it is a local one stipulated, for 
instance, by the passage of some superpowerful torrent 
of cosmic particles or, finally, whether it is a result of a 
combination of both alternatives. 

We may hope that a further detailed statistical inves- 
tigation of plates will bring some answer to these questions. 

1L. Myssowsky and P. Tschishow, Zeits. f. Physik 44, 408 (1927); 
A. Jdanov, J. de phys. et rad. 6, 233 (1935). 

2See M. Shapiro, Rev. Mod. Phys. 13, 58 (1941) for a summary. 

3A, Jdanov, Bull. Acad. Sci. U.S.S.R. 4, 266 (1940) ; Comptes rendus 
Acad. Sci. U.S.S.R. 28, 109 (1940). 

+A. Filippov, A. Jdanov, and I. Gurevich, Comptes rendus Acad. Sci. 


U.S.S.R. 18, 181 (1938); J. Phys. U.S.S.R. 1, 51 (1939); A. Jdanov, 
Comptes rendus Acad. Sci. U.S.S.R. 20, 645 (1938); 22, 163 (1939). 


6 If one takes into consideration only those disintegrations and single 
tracks where the range is more than 10 cm of air equivalent. 
A. Jdanov, Nature 143, 682 (1939). 








The Method of Shower Anticoincidences for 
Measuring the Meson Component of 
Cosmic Radiation 


VIKRAM SARABHAI 


Department of Physics, The Indian Institute of Science, 
Bangalore, India 


January 24, 1944 


HE meson intensity measured in counter experiments 

where absorbers are used to filter out the soft elec- 
tronic component refers only to fast mesons. The slow 
mesons which get cut out by the absorber can so far be 
estimated only indirectly, as, for instance, by the method 
adopted by Auger! and Greisen.* There is evidence, how- 
ever, that the proportion of slow mesons increases rapidly 
with altitude, and it would be of great advantage if there 
were a direct method for measuring the total meson 
intensity. Bhabha* has recently suggested on grounds of 
reducing the weight of lead carried by balloons measuring 
the meson intensity in the stratosphere that the production 
of secondaries by the electronic component can be used 
to advantage in cutting them out. But this principle could 
be used for the more important purpose of bringing for the 
first time the slow mesons in the field of direct experi- 
mental observations. While a great proportion of cosmic- 
ray electrons in the atmosphere are already associated with 
other ionizing particles, the number so associated can be 
increased by a lead plate of thickness corresponding to the 
maximum of the Rossi curve. By means of an anticoin- 


cidence arrangement working along with a vertical counter © 


telescope, it is possible to exclude events in which associated 
particles arrive in the telescope. This can be done with side 
anticounters arranged as shown in Fig. 1a, but here the 
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Fic. 1. a. Arran ent showing anticounters outside the main cone 
of measured radiation. b. Arrangement used for measuring shower anti- 
coincidences. c. Basic circuit for registering shower anticoincidences. 



































main cone of measured radiation is not utilized for detect- 
ing the showers and those that are generated in the lead 
above on account of their narrow angular spread might 
avoid detection. The actual arrangement used is shown in 
Fig. 1b, where the sets of counter II and III are used to 
register the showers. A shower coincidence II, III gives 
rise to an antipulse which is fed to the coincidences of the 
main telescope formed by counters I, IV. The lead placed 
in position X is to increase the associated soft particles, 
and the lead placed in position Y is to study the absorption 
of the radiation without interfering with the efficiency of 
shower anticoincidences. The basic circuit used has been 
shown in Fig. 1c. 

Preliminary results obtained with this arrangement of 
the total intensity given by coincidences I, IV and the 
meson intensity given by shower anticoincidences I, IV- 
(II, III) at various altitudes up to 13,900 ft. in Kashmere 
illustrate the working of the arrangement. Figure 2 shows 
the absorption curve at two altitudes of these two counting 
rates with lead placed at positions X and Y. The difference 
between the two curves diminishes rapidly with increasing 
thickness of lead absorber as the electronic component gets 
progressively absorbed; and while the slope of the total 
intensity curve increases with diminishing thickness of 
lead, the slope of the meson intensity curve varies in an 
opposite direction. The tendency of the shower anticoin- 
cidence curve to become horizontal for zero thickness of 
lead is a striking indication of the fact that it measures 
chiefly the meson intensity. The altitude intensity curves 
from readings at 4 elevations in Kashmere at geomagnetic 
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Fic. 2. Curves showing for two altitudes the absorption of the otal 
intensity measured by coincidence rate I, IV (full line) and the meso: 
intensity measured by shower anticoincidence rate I, IV-(II, Il) 
(broken lines). 


latitude 25°N show that above 10,000 ft. the meson 
intensity increases less rapidly than the total intensity and 
the atmospheric shower intensity. At 13,900 ft. the pro- 
portion of mesons was 64 percent of which about 18 percent 
consisted of slow mesons, but the accuracy of these 
estimates is not high. A full report of the results obtained 
at Kashmere and repeated with the same arrangement at 
Bangalore will appear shortly.‘ 

The errors in this experiment might be due to three 
causes. Firstly, due to the contribution of side showers to 
the double coincidences that were measured. This would 
give an under-estimate of the proportion of electrons. 
Secondly, due to electrons of low energy which do not 
produce large enough showers to be detected and cut out 
by the anticoincidence arrangement, and lastly, due to 
mesons which are associated by other mesons or knock-on 
electrons and which get cut out by the anticoincidence 
arrangement. These two errors will tend to give an over- 
estimate of the proportion of slow mesons. With the help 
of an accurate experiment that is now being performed with 
quadruple coincidences measuring the vertical intensity, 
it is hoped to correct for these errors. 

I am indebted to Professor Sir C. V. Raman for supprt 
and assistance, and to Professor Bhabha for helpful dis- 
cussions. The Kashmere trip would not have been possible 
but for the willing help given by the state officials and 
friends. 

1P. Auger, Phys. Rev. 61, 684 (1942). 
2K. Greisen, Phys. Rev. 63, = (1943). 


* Bhabha, Proc. ind. Acad. Sci. (A), in press. 
4 Sarabhai, Proc. Ind. Acad. Sci. (A), in press. 





Configurations of a Binary Mixed Crystal 
, R. EIsENSCHITZ 


Davy Faraday Laboratory of the Royal Institution, London, England 
January 31, 1944 


N a recent paper Ashkin and Lamb! applied the matrix 
method of statistical mechanics to mixed crystals, 
especially to the order-disorder equilibrium. In their theory 
any configuration of the crystal is specified as usually by 


the kind of atom (A or B) by which every lattice point js 
occupied. The matrices by which the correlations in the 
lattice are determined appear to be “degenerated” at low 
temperatures; the quotation mark indicates that it is the 
absolute and not the algebraic value of the highest proper 
value which matters. By this method information is ob. 
tained on the superstructure, whereas the evaluation of 
the specific heat curve seems to be difficult. 

In this note an alternative specification of configurations 
is discussed which is equivalent to the above specification, 
In a simple cubic lattice of a binary alloy such as B-brass, 
there are two possible values for the bond energy ¢ 
namely, e=0 and e=g¢ corresponding to bonds between 
atoms of different kinds and the same kind, respectively, 
¢ is assumed to be positive. The configurations may ae. 
cordingly be specified by the pattern of bonds, i.e., by the 
value of ¢ for every bond in the lattice. Not every pattern 
of bonds corresponds to an atomic configuration; only 
those patterns are admissible in which there is an eveg 
number of “‘O bonds” on every closed contour of bonds jp 
the lattice. Every pattern which complies with that con. 
dition corresponds to two atomic configurations which are 
transformed into each other by interchanging all A atoms 
with B atoms and vice versa. The set of all admissible 
patterns corresponds to the set of all atomic configurations 
in which not only the positions of A and B atoms but also 
their ratio is varied. 

Instead of making statistics over atomic configura 
tions, it may be made over the admissible patterns of 
bonds. The central problem is the enumeration of those 
admissible patterns for which the number of “¢ bonds” 
and accordingly the configurational energy is given. Ina 
square of four nearest neighbors those patterns are admis- 
sible in which the number of O bonds is 0 or 2 or 4. There 
are altogether 8 patterns of this kind. If a pattern of bonds 
in the lattice is admissible, then every square necesssarily 
exhibits one of these 8 patterns. This necessary condition 
is also sufficient. Assuming that the number of O bonds 
is even in every square of four nearest neighbors, it is 
proved in the following that this number is even on every 
closed contour of bonds in the lattice. 

In proving this, at first a square net in the plane is con- 
sidered. Set g= /(1—«/g)ds a contour integral taken 
counterclockwise over the circumference of one square. 
This integral is equal to a sum of four terms which are 0 
for every ¢ bond and +1 for every O bond. Let Q be the 
sum of all g taken over an area which is enclosed by any 
arbitrary contour of bonds in the net. Then all bonds in 
the interior of this area contribute 0 to the sum of contour 
integrals. Q is equal to the contour integral /(1—«/¢)ds, 
which is taken over the boundary and equal to a sum of as 
many terms as there are bonds in the boundary. Since, 
according to premises, every square contributes in even 
number to Q, Q is even and so is the sum taken over the 
boundary. As that sum consists only of terms equal to |, 


‘0, —1, the number of its non-vanishing terms is neces 


sarily even. There is accordingly an even number d 
O bonds on the boundary, whatever the contour may be 
like. 
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In a lattice of three dimensions an arbitrary closed 
contour of bonds is considered to be the boundary of a 
large number of surfaces which are formed by squares 
between nearest neighbors in the lattice. Two of these 
surfaces define a closed surface and accordingly the direc- 
tion of an outer normal for each square on either surface. 
Having assigned on one surface an outer normal to each 
square, ¢ is defined as contour integral taken counter- 
clockwise over one square as seen from the outer normal 
and Q is defined as the sum of all g taken over the surface. 
According to the above argument, the number of O bonds 
on the contour is then seen to be even. 

So far the number » of A atoms in the lattice of m atoms 
is regarded as variable so that the configurational partition 
function is 


P= piv), 
v=0 


where p(v) is the partial sum over the configurations with 
» A atoms. The partition function of an equimolecular 


alloy is 


P*= 2 p(v)-exp (— ap), 
pad 


where the parameter a is determined by 
n/2=—d log P*/da. 


Since p(v) =p(n—v), it follows from the former equation 
that a=0 for all temperatures. The specific heat curve of 
an equimolecular alloy is accordingly equal to the specific 
heat curve of an alloy with variable composition. 

The writer* recently applied the matrix method to the 
statistical calculation of the specific heat. The con- 
figurations of the lattice are specified in this calculation 
by the pattern of bonds and only such patterns are 
admitted in which every square has an even number of 
O bonds. It appears that in this case the matrices show 
no “degeneracy.”” While no information is obtained on the 
superstructure, the specific heat is evaluated for all tem- 
peratures; it has no singularity. A note on the correspond- 
ing result for a two-dimensional net has been published 
previously.* 

1 J. Ashkin and W. E. Lamb, Phys. Rev. 64, 159 (1943). 


?In print in the Proceedings of the Royal Society. 
3R. Eisenschitz, Nature 147, 779 (1941). 





Relativity and Supplementary Indeterminacy 


G. WATAGHIN 
Department of Physics, Sio Paulo University, Sio Paulo, Brazil 
February 5, 1944 


N recent years experimental evidence has been ac- 

cumulated concerning the production of groups of 
mesotrons in one collision and processes in which the 
hypothetical neutrinos take part.! It was pointed out? that 
in view of the irreversibility of many of the processes ob- 
served with cosmic rays the statistical mechanics of an 
assemblage of particles having an average energy per par- 
ticle E=200 mc? is quite different from the usual one. For 
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instance, the introduction of temperature and entropy is 
not allowed in the cases in which statistical equilibrium is 
impossible. New limitations arise because of the following 
circumstances: In a collision accompanied by the produc- 
tion of several mesotrons, the rest mass of the group 
changes, and the control of the conservation laws of mo- 
menta and spins appears limited by the decay processes and 
fluctuation of the number of created mesotrons and 
neutrinos. The wave properties of particles, based on dif- 
fraction and interference phenomena, are also altered and 
probably unobservable in high energy collisions, in ac- 
cordance with the non-linear character of the laws (differ- 
ential equations, interaction operators) governing such 
collisions. The representation of high frequency photon- 
field (for hy >200 mc?) as a system of harmonic oscillators, 
and the statistical laws of Planck, Fermi, Bose at extremely 
high temperatures are invalidated for similar reasons. From 
the foregoing considerations, it follows that the behavior of 
high energy particles (observed as cosmic rays in a terres- 
trial frame of reference) is qualitatively different from the 
usual one, in accordance with the idea of supplementary 
indeterminacy .* 

In the theory of relativity we must distinguish the postu- 
late of general covariancy referring to an arbitrary trans- 
formation of coordinates (mathematical symbols used to 
denominate the events in 5S,) from the equivalence of 
different physical reference frames, each individualized by 
a system of measuring devices (diffraction gratings, micro- 
scopes). Obviously our possibilities of building a reference 
frame corresponding to a given mathematical transforma- 
tion are limited (we cannot use rigid bodies or apparatus 
moving with v>c). Let us consider a terrestrial frame of 
reference (neglecting the rotation) as a local “‘co-moving”’ 
frame in the Friedman universe in expansion. The ve- 
locities of the receding nebulae and the cosmic rays are 
distributed around the earth with spherical symmetry. Let 
us consider another Lorentz observer moving with an ultra- 
relativistic velocity 8c [(1—6*)-4>1000]. Any particle of 
the apparatus used by this second observer will constitute 
for the terrestrial observer a high energy particle capable of 
producing a cascade shower or a mesotron shower. Now we 
introduce the assumption that the “ultra-relativistic’’ 
frames of reference are not equivalent to the co-moving 
ones because the apparatus moving with the velocity of 
high energy cosmic rays could not function or could not be 
used as measuring devices in the actual world. Thus we can 
exclude from our relativity postulate such non-allowed 
reference systems. 

In this way we overcome the incompatibility between the 
relativity and the theory of a world in which a universal 
length marks the lower limit for measurement of intervals 
and a supplementary indeterminacy exists in the region of 
high energies and momenta. 

1See cloud-chamber photographs of groups of three, five, seven, 
mesotrons; experiments with counters on mesotron showers produced 
in air, lead, and also in clay, underground; high altitude researches on 
the production of protons, neutrons, and mesotrons in lead and paraffin. 
Some of the secondary protons and neutrons behave !ike primaries and 
thus generate mesotron showers. 

2G. Wataghin, Phys. Rev. 63, 137 (1943), and 64, 248 (1943). 


? Nature 142, 393 (1938). 
4 Zeits. f. Physik 66, 650 (1930); Comptes rendus 207, 421 (1938). 
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On the Spin of the Meson 


S. K. CHAKRABARTY AND R. C. MAJUMDAR 
University College of Science and Technology, Calcutta, India 
February 14, 1944 


HE phenomena of large bursts under thick layers of a 
material consisting of more than 200 particles, which 
are produced almost entirely by mesons of energies greater 
than 10" ev provide a very feasible means for the determi- 
nation of the spin of the meson. Christy and Kusaka! calcu- 
lated the frequency of large bursts produced by mesons and 
concluded that the results for a meson of mass 177 times 
that of electron and spin zero are in satisfactory agreement 
with the observation of Schein and Gill.* Similar calcula- 
tions made by one of us (S.K.C.)? showed that the con- 
clusions of Christy and Kusaka cannot be maintained and 
that the observational results fitted well with the results of 
calculations for a meson of spin 1 when one took the 
bremsstrahlung cross section, after taking into considera- 
tion the effect of radiation damping in it, as given by 
Wilson.‘ The reasons for arriving at such opposite con- 
clusions were discussed in (C) (see reference 3). It now 
appears to us that the analyses of Wilson contain some 
errors, and the correct expression for it may somewhat 
modify the results obtained (C) in this particular case. By 
introducing properly the effect of the radiation damping on 
the emission of photons by mesons, we obtain the differ- 
ential cross section ¢(¢)de for the emission of a photon of 
energy lying between eW and (e+de) W by a meson of spin 1 
and energy W, as given by the following expressions: 


o(e)de= Gde[ Cofo(a)—Cifi(a)+L(W, 6)], for a<i (1) 


and 
(e)de= Gael Co F log 2+fo(*) } 


— C742 logta—fi(+) 4-200, 9] or 431 @) 


where 
L(W, b) =2C, log? (5) — Cofo(b) + Cifi(b) 


(3C2+8)g8 


3b +2C, log (1+5?) 


— Co tan of log (5) — BEr¥ 2008 26)8 


fola)=a(1-F+5- - ‘), 


+ 


fila) =a*(1-5+5,— .* ‘), a=WB/«MC?, 


2—2e+72 


C= 12B 


—(C2+4)8, 
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_ (34—34e+ 7e)e 


ea, ~~ * 
_ (14—14e+2)e . 
C= 24(1—«)? ”. 


_& 


= bo = ‘ie = . a 
B= (5/288)'a~9.63X10", b= 57>, 


G=aZ*7F(m/M)*, a=e/hC. 


«x is 1 or 5Z!/m according to whether we take for the lowest 
value of the impact parameter the Compton wave-length 
(Wilson) or the nuclear radius (Christy and Kusaka), The 
cross section for the case of spin 0 will, however, not be 
appreciably affected by the inclusion of the effect of radia. 
tion damping. Whether the frequency of bursts calculated 
with (1) and (2) will still compare favorably’ with that 
obtained in (C) cannot be stated without making a long and 
tedious calculation. 

While this calculation is in progress we meanwhile point 
out that the frequency of bursts produced by the pseudo. 
scalar mesons of spin 0 can be made to agree with the 
observation of Schein and Gill by a suitable adjustment of 
its mass. From Eq. (18a) in (C) the total probability B(N) 
of getting a burst containing more than N particles below 
an infinitely thick layer of lead, for different values of », 
the ratio between the mass of the vector and pseudoscalar 
mesons, have been calculated. The mass of the vector 
meson has been taken as in (C) to be 177 times the electron 
mass. The results obtained for two different values of y are 
given in Table I. On comparing these values with the 


TABLE I. Values of B(N)X10° for a pseudoscalar meson. 











N 
™ 100 200 


* ™\ 
2.5 355.5 99.75 28.00 13.41 7.976 5.191 3.733 2.203 
3.0 527.0 147.4 41.30 19.76 11.74 7.639 5492 3.23% 
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observations of Schein and Gill we find that the agreement 
is satisfactory when y=~2.5. Moreover, our calculations 
show that in order to obtain the correct sign and magnitude 
of the quadrupole moment of the deuteron in the theory of 
Mgller-Rosenfeld and Schwinger,* we shall have to take 
y ~2.5. (Pauli and Kusaka’ find this ratio to be greater 
than 2.) These results therefore speak in favor of the theory 
of Mgller and Rosenfeld as modified by Schwinger. The 
analysis will be shortly published in detail. 

1 R. F. Christy and S. Kusaka, Phys. Rev. 59, 405, 414 (1941). 

2M. Schein and P. S. Gill, Rev. Mod. Phys. 11, 267 (1939). 

3S. K. Chakrabarty, Ind. J. Phys. 16, 377 (1942). To be quoted as (C), 

4A. H. Wilson, Proc. Camb. Phil. Soc. 37, 301 (1941). 

5C, Moller and L. Rosenfeld, Nature 144, 476 (1939). 


6 J. Schwinger, Phys. Rev. 61, 387A (1942). 
7 W. Pauli and S. Kusaka, Phys. Rev. 63, 400 (1943). 
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